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Preface 


These exercises have been developed to fit the content of the text book Drilling Fluid Engineering at 
www.bookboone.dk . The understanding of the physics and mathematics of the processes has been in 
focus of both the textbook and the exercises book. Many practical applications have also been created and 
entered into the collection of exercises. Most of the exercises have been solved and corrected by students 
in the corresponding course at the Department of Petroleum Engineering and Applied Geophysics at 
NTNU in Trondheim. If the readers have any comments that could improve the exercises, please contact 
me at pal.skalle@ntnu.no . Any such comments will be worked into the next years issue of this book. 

Pal Skalle 

Trondheim, oktober 2015 
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1 Fluid Properties 

We moved exercises 1.1-1.4 to other chapters. 

1.5 Rheology control 

a) Will YP, PV and p eff be influenced by the addition of barite? 

b) Why is lye (NaOH) added to the drilling fluid? 

c) Define polymers and the purpose of adding them to the drilling fluid. 

d) Define pseudo plastic, thixotropic and rheopectic fluid behaviour. 

e) Why does viscosity of water increase when Bentonite is added? 

1.6 Rheology control 

a) Explain how dispersed Bentonite is able to contain up to 18 times its own volume of distilled 
water. Why is it that the water-holding effect will be reduced when salt is added to the water? 

b) Explain the reason behind the non-Newtonian behavior of Bentonite suspensions. 

Out on a drilling rig the questions asked are of practical nature: In the upper wellbore section seawater 
is often used as drilling fluid. If the viscosifying effect drilled-through clay does not produce the proper 
viscosity, the addition of Bentonite to the water has to be considered. Use Figure 1-6 and assuming that 
the quality of the drilled-out clay corresponds to Premium Drilling Clay. Assume that ROP is 35 m/hr 
while using a 26” bit. The pump rate is 3 000 1/min. The required mud viscosity must at least be 15 cP. 

c) Will the formation provide the required viscosity? 

d) What is the yield (m 3 of mud/ton solids) of Wyoming Bentonite? (p Bentonite = 2.4 kg /l) 

e) What is the water density (originally 1.0 kg/1) after addition of Bentonite, when an effective 
viscosity of 50 cP is the upper boundary? 
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Figure 1 -6: The ability of different solids to produce viscosity. 

1.7 Flocculation 

Bentonite and polymers are the dominating viscosity agents. Bentonite is still widely used and we 
have to understand its behavior properly. Bentonite behaves different from other additives; it swells 
and flocculates. 

a) Why does Bentonite flocculate (weak flocculation)? 

b) What will happen if untreated water based mud is used while drilling through the cement in 
the casing shoe area? Explain what happens to the mud (strong flocculation) and the respective 
operational consequences. Sketch the flow curve of the mud before and after having drilled 
through the cement. 

c) Name 3 factors which enhance flocculation of Bentonite. 

d) Why does WBM behave shear thinning? 

e) Why does WBM behave thixotropic? 

f) Explain the difference in flocculation tendency of pure edge-to-face and of cross-linking caused 
by external agents. 
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Fluid Properties 


1.8 Mud contamination 

a) What is a contaminated mud, and how are the drilling fluid parameters restored? 

b) The geologist expected that layers of silty anhydrite would be penetrated at a vertical depth 
of 1600 m. The mud engineer was therefore told to make measurements every 15 min. of the 
returning mud as shown in the table below. After penetrating the anhydrite he observed that 
the viscosity of the mud, a dispersed WBM system, started to rise and became abnormally 
thick. Drilling continued and after a few hours of drilling/pumping, the viscosity fell back to 
a lower level than the original viscosity. 


Parameter 

Shear stress at RPM of 
600 300 

100 

P 

V filtrate 

Unit 



(lb/100ft 2 ) 


kg/l 

ml / 30 min 


0900 

42 

28 

16 

1.31 

7 

[ Time 

0915 

41 

27 

17 

1.31 

6 

0930 

68 

54 

37 

1.31 

18 


1145 

31 

17 

5 

1.30 

7 


Together with the mud engineer you are responsible for the maintenance of the drilling fluid 
program. Explain changes in the recorded parameters observed at 0930 and at 1145. Suggest 
countermeasures against these changes. 

c) Assume two clay suspensions are flocculating for two different reasons during drilling 
operations; 1. Edge-to-face. 2. Calcium attack. How do the two clay suspensions behave 
Theologically? Make a sketch of the in-situ shear stress vs. time; during drilling into the Ca ++ 
containing layer. Assume you drill in the contaminated zone for 10 min. Then turn off the 
pump for 10 min. Continue drilling for 10 min till you are out of the contaminated zone. 

1.9 Flocculation 

a) Find necessary YP to keep a spherical particle suspended in a mud of = 1.1 kg/1. The 
particle has these characteristics: 

d p = 5 mm 
p p = 2.3 kg/1 

Similarly, find out what is the maximum size the particles can be kept suspended at when 
YP = 15 Pa. 

b) At 12:00 the ROP became very low and it was decided to change the bit. A 15 min. stop in 
the operation was made before tripping-out from 2 100 mMD was initiated. When the first 
pipe was broken mud spilled out on the drill floor. Could this spill have been prevented? How 
high up can the string be hoisted before gravity pulls the mud down in the following situation: 

x gd 15 min. = 30 lb/100 ft 2 (14.4 Pa) 

Inner diameter of the drill pipe = 4.127" (104.85 mm) 
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1.10 Fluid additives 

a) Define the different concepts and explain their relevance for drilling fluids (e.g.; viscosifyer). 


Anhydrite 

Lignosulphonate 

Caustic soda 

Lignite 

CEC 

MBT 

Chalk 

PAC 

CMC 

Pre hydrated 

Colloid 

PHPA 

Dispergator 

SAPP 

Deflocculators 

Sodium Sulphate 

Gypsum 

Starch 

HEC 

Xanthan 

NaOH 



b) How does the additive called drag reducer reduce turbulent pressure so dramatically? 

c) What significance does the K + concentration have for the shale? 

d) When drilling into swelling clay and swelling shale, problems like sloughing shale and stuck 
pipe may occur. Explain what happens to the mud. How and why do you convert Bentonite 
mud into gyp mud? 

1.11 Fluid additives 

a) How many moles/liter of hydroxyl (OH) concentration is required to change the pH of a 
drilling fluid from 7.5 to 11? 

b) How much caustic soda (weight per liter) will be required to increase the pH in question a)? 

c) Why potassium hydroxide (KOH) is often preferred to sodium hydroxide (NaOH) in controlling 
the pH of mud? 

How does drilling fluid achieves the following functions: 

d) Lift cuttings from the bottom to the surface 

e) Releases cuttings at the surface 

f) Cools and lubricates the drill bit and the drill stem 

g) Prevent blowouts 
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Rheological models 


2 Rheological models 

To simplify the evaluation of drilling fluids out in the field, the simplified Bingham field method was 
developed. In this collection of exercises we distinguish between the simplified field and the standard 
method of determining rheological model constants. 

The simplified field method is applicable in conjunction with the Fann viscometer. 

For the standard method the Fann readings are converted to SI units and multiplied with the factor 
1.067 (see SPEs Applied Drilling Engineering textbook Appendix A, eqn. A-6b). 

Note also that conversion factors are presented in Chapter 6 in present book. 




MAERSK 


I joined MITAS because 
I wanted real responsibility 


The Graduate Programme 
for Engineers and Geoscientists 

www.discovermitas.com 


Real work 
International opportunities 
Three work placements 


0 


a 


I was a construction 
supervisor in 
the North Sea 
advising and 
helping foremen 
solve problems 


14 

Download free eBooks at bookboon.com 









Exercises in Drilling Fluid Engineering 


Rheological models 


2.1 Bingham/Power law. 

In the laboratory the data in Table 2-1 were obtained (0 is the dial reading in the viscometer): 


Parameter 

RPM 

y 

0 

T 

T 

Unit 


(s- 1 ) 

(-) 

(lb/100ft 2 ) 

(Pa) 

Data 

600 

1022 

106.4 

112.8 

54.0 

300 

511 

75.0 

79.5 

38.1 

100 

170 

42.4 

44.9 

21.5 

6 

10 

10.0 

10.6 

5.1 

3 

5 

6 

6.3 

3.0 


Table 2-1: Rheological data. 


The same data are presented as a flow curve in Figure 2-2. 



Figure 2-1 : Flow curve of the rheological data from the table 2-1. 

a) Find rheological constants for the two rheological models Bingham and Power-law. For 
Bingham model use both field and standard method. 

b) Which of the two models in question a) fit the shear stress best at 100 RPM. 

c) Which of the three models, Newtonian, Bingham or Power-law, would give the best answer 
on basis of the given rheology while pumping 1000 1/min through a pipe of 10 cm ID. Hint: 
Check theoretical vs. measured shear stress. 
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2.2 Bingham/Power-law 

Rheological data are tabulated and presented graphically in Figure 2-2. Note that shear stress, t, has 
been multiplied by 1.06 before converting readings, 0, to Si-units. 




Figure 2-2: Graphical representation of rheological data (flow curve). 


10 100 1000 10000 
Shear Rate (s ' 1 ) 



a) Select the best 2-data-point-rheology model, either Bingham or Power-law (at an viscometer 
speed of 100 RPM). Verify selection. 

b) Observe the log-log plot. It represents a typical clay-dispersed system. 

Why do the two data points of the lowest shear rate deviates from the straight line made on 
basis of the upper four data points? 

2.3 Bingham/Power-law. Regression 

After measuring the rheology of the fluid it is always useful to plot its flow curve. The following data 
are obtained: 


Text 

Speed 

Reading 

Shear stress 

Symbol 

Y 

e 

T 

Unit 

(rpm) 

(s 4 ) 

(-) 

(lbf/100 ft 2 ) 

(Pa) 


600 

1022 

60.4 

64 

30.6 


300 

511 

39.6 

42 

20.2 

Data 

200 

340 

32.1 

34 

16.3 


100 

170 

24.5 

26 

12.5 


6 

10 

14.2 

15 

7.2 


3 

5 

9.4 

10 

4.8 
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Rheological models 


a) Plot x vs. y for three rheological models Newtonian, Bingham, Power-law according to 

Field procedure (2 data points) (only for Bingham) 

Standard procedure (2 data points) (for all models) 

Regression (6 data points). Use Excel linear regression. 

Find the constants of the three models. Determine which of the models are best fitted to the 
readings at high shear rates, i.e. for determining pressure loss in nozzles and inside the drill 
string (300 rpm and higher), and what model fits best to the lower shear rates, valid for annulus 
or flow (100 rpm and lower). 

b) An exercise without calculations: Will plug flow occur in this mud system? 

2.4 Effective viscosity 

a) Rheology data: 


Rotational speed 

Dial Reading 

600 rpm 

43 

300 rpm 

30 

200 rpm 

23 

100 rpm 

16 

6 rpm 

8 

3 rpm 

7 




Because achieving your dreams is your greatest challenge. IE Business School’s Master in Management taught in English, 
Spanish or bilingually, trains young high performance professionals at the beginning of their career through an innovative 
and stimulating program that will help them reach their full potential. 

Choose your area of specialization. 

Customize your master through the different options offered. 

Global Immersion Weeks in locations such as London, Silicon Valley or Shanghai. 

Because you change , we change with you . 


www.ie.edu/master-management mim.admissions@ie.edu f # In YnTube ii 


Master in Management 


17 

Download free eBooks at bookboon.com 




















Exercises in Drilling Fluid Engineering 
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A useful exercise is to plot the effective viscosity (apparent Newtonian viscosity) as a function 
of shear rate. The non-Newtonian and the shear thinning effect will then appear clearly 

b) The mud is flowing in a 1000 m long pipe with inner diameter of 4 in, flow rate is 6 000 1/ 
min (two pumps) and density is 1.1 kg/1. Use the Power-Law to estimate p eff at this flow rate. 

c) Determine the pressure loss in the pipe (Power-law). 

d) Why are drilling fluids often so well suited to the Bingham model? 

2.5 All models 

The flow rate is 2 500 1pm in a 1000 m long pipe with an inner diameter of 10 cm. Rheological data 
points are given below. The mud density is 1.1 kg/1. 


Shear rate (s' 1 ) 

Shear stress (Pa) 

1022 

55 

511 

40 

340 

35 

5 

19 


a) Find the rheological constant just for Bingham and Power-law (use only upper 2 data points). 

b) For Herschel-Bulkley model, discuss three different ways of obtaining the constants without 
calculation. 

c) Show that the effective viscosity for Bingham fluids is: 

r 0 d 

Veff = Vpl + ~^— 

d) Explain why the Bingham field model is so useful for evaluating mud behavior. 

2.6 All models. Regression 


Symbol 

Y 

RPM 

0 

T 

T 

Unit 

Hz 

rpm 

- 

lb/100 ft 2 

Pa 


1022 

600 

50 

53 

25.4 


511 

300 

34 

36 

17.2 


340 

200 

26 

27.6 

13.2 

Data 

170 

100 

17 

18 

8.6 


102 

60 

13 

13.8 

6.6 


51 

30 

9 

9.5 

4.5 


10 

6 

4 

4.2 

2.0 


5 

3 

3 

3.2 

1.5 


a) Use the Fann-viscometer data above to determine mo del-constants for the first 2 models listed 
below through standard (2 data points) procedure. 
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Rheological models 


b) Perform linear regression procedure (use Excel spread sheet). Make a plot of the two first 
models for 2, 4, 6 and 8 data points. Linear regression is in fact possible only for those two 
models. The lower three models are presented just to give an overview of rheological models. 


Bingham 
Power law 

Herschel Bulkley (H-B) 
Collin-Graves (C-G) 
Robertson-Stiff (R-S) 
Casson 


x = x + ii *y 
x = K* y n 
x = x + K y n 
x = (t q + K y n ) (l-e _,3y ) 
x = K(y o + y ) n 

r= yr 0 + vV’yJ 


c) To solve H-B, use two methods: 
cl) Elimination and iteration. 
c2) Non-linear regression. 


The latter procedure is presented in Chapter 10.6 in the Drilling Fluid Engineering Text book. 


2.7 All models 


Parameter 

Speed 

Shear Stress 

Unit 

Hz 

lb/100 fl 2 

Pa 


1022 

52.2 

25 

Data 

511 

35.5 

17 


340 

25.1 

12 


170 

13.6 

6.5 


10 

6.3 

3.0 


5 

4.2 

2.0 


Which of the following three models, Power-law, Bingham and Herschel-Bulkley would you select for 
estimating x at y = 10 Hz. Apply the two and three upper data points. Use Field and standard procedure 
for the Bingham 2-data point model. For the Herschel-Bulkley model, use the x 5 reading as x . K and n 
are found from Power-law. 


r = r 0 +//-y 

T = K-/ n 

r = T 0 +K-r n 
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Drilling fluid dynamics 


3 Drilling fluid dynamics 

3.1 Velocity profile. Continuity equation 

The incompressible steady state flow between two parallel plates with breath b in Figure 3.1 is initially 
uniform at the entrance; v = v = 8 cm/s. Downstream the flow develops into the parabolic laminar profile 
v(z) = az (z A - z), where a is constant and z A the plate distance. If z = 4 cm, what is the value of v ? 

v/v 0 7 0 r 0 max 



Figure 3-1 : Flow data. 


3.2 Velocity profile. Momentum flux 

a) The fully developed laminar pipe-flow velocity profile is expressed as: v z (r) = v max (l - r 2 / R 2 ), 
v e = 0 , v r = 0. z indicates here the axial direction: This is an exact solution to the cylindrical 
Navier-Stoke equation. Neglect gravity and compute the pressure distribution in the pipe; 
p(r,z), and the shear-stress distribution; x(r,z), using R, v max and p as parameters. Why does 
the maximum shear occur at the wall? 

b) For flow between parallel plates, compute bl) wall shear stress and b2) the average velocity. 
From the Text book, Chapter 4, we find that v z (y) = - dp/dz h 2 /2 ju (1 - y 2 /h 2 ). The most 
general definition of shear stress is given by: r = t z = p -I- 



Figure 3-2: The geometry of pipe flow in the z-direction. 6 = 0 in this exercise. 


3.3 Velocity profile 

a) Discuss the meaning of this expression, its assumptions etc. 


dp Id dv 

— =-( r ju —-) 

dx r dr dr 
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Drilling fluid dynamics 


b) For stationary non rotational and laminar flow of Newtonian fluids in circular horizontal pipes 

(z = x), show that v(r) = dpidx ^ _ r ^ 

A fd 

c) Determine average velocity. The absolute velocity is largest in the center of the pipe. Max velocity 
compared to the average velocity is forming an expression of the axial dispersion when one 
fluid is displaced by another. Find this expression. 

d) Find wall shear stress and average pipe velocity when the pressure loss is recorded to be 0.9 
bar along a 1 000 m long pipe. Radius is 5 cm and viscosity 53.8 cR 

3.4a Pressure loss vs. rheology 

Water, assumed incompressible, flows steadily through a pipe of constant diameter 2R. The entrance 
velocity is constant, u = u o , and the exit velocity approximates turbulent flow, u = u max (1 - r/R) 1/4 . Fluid 
viscosity is p. Determine the average velocity and the shear stress at the wall during turbulent flow. 
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Drilling fluid dynamics 


3.4b Pressure loss vs. rheology 

Use the viscometer readings from exercise 2.6. 


a) Determine the pressure loss pr. 1000 m in a 10 cm ID pipe at a flow rate of 1000 1/min. and a 
fluid density of 1000 kg/m 3 . Use three rheological models. The observed pressure loss at these 
circumstances were 4.5 bar. 

b) Calculate shear rates in the pipe for Bingham, Power-law and Newtonian model. Read shear 
stress from the flow curve, and determine pressure drop through the universal pressure 
loss model. 

L 

Ap = Tw ’ d ’ 4 

l 

c) Show that wall shear stress can be expressed as T w = - 
for laminar, annular pipe flow, and that fi e ff = fi p i + 

d) What effect has entrance length of a uniform pipe on estimated pressure loss? 

3.5 Pressure loss vs. rheology 

Mud is pumped at a rate of 800 1/min with these rheological data: 


n ■ 

M 


T 'dhydr 


8V 


RPM 

0 

x fib/100 ft 2 l 

600 

66.0 

70 

300 

47.2 

50 

100 

25.5 

27 

6 

9.4 

10 


a) Find Ap p . pe in a 1000 m long pipe of d = 0.109 m for a field-Bingham fluid. Fluid density is 
1100 kg/m 3 . 

b) Which rheological model, Newtonian or field-Bingham, is better suited for pressure loss 
estimation when the actual Ap p . pe was recorded to be 0.7 MPa. Apply the universal pressure 
loss model; Ap = 4 x L / d. 
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3.6 Pressure loss. Power-law 

a) Derive the laminar pressure loss expression in pipes for Power Law fluids from the force balance. 

b) Show that: 

pv_d = N 

R e generalized 

when p eff for a Power-law fluid is applied. 

c) Show that the Reynolds number for a Power-law fluid increases as the inner pipe diameter of 
the annulus decreases, while for a Newtonian it decreases. Apply annular diameter in terms 
of d, , = d - d., and let flow rate be constant. 

hydr o 1 

3.7 Pressure loss. Turbulent. Energy equation 

Oil of the density p = 900 kg/m 3 and a kinematic viscosity v = 0.00001 m 2 /s, flows at a rate of 0.2 m 3 /s 
through a 500 m new cast-iron pipe with a diameter of 200 mm and a roughness of 0.26 mm. Determine 
the head loss. 


3.8 Pressure loss vs. flow rate 


Make a graph of pressure loss vs. flow rate in a 1000 m long pipe with inner diameter of 10 cm. Rheological 
data points are given in Exercise 2.6, and the mud density is 1.1 kg/1. Select the Power law model. 



AACSB 


ACCREDITED 


Excellent Economics and Business programmes 

‘u&r 


university of 
groningen 




www.rug.nl/feb/education 


“The perfect start 
of a successful, 
international career.” 

CLICK HERE 

to discover why both socially 
and academically the University 
of Groningen is one of the best 

places for a student to be 


23 



Download free eBooks at bookboon.com 



Exercises in Drilling Fluid Engineering 


Drilling fluid dynamics 


3.9 Pressure loss. Use field data to evaluate model 

Prior to a pre flush/cementing operation the driller performed pressure tests to verify theoretical 
pressure estimations. Previously comparisons between the actual pressure readings during drilling with 
theoretically estimated pressure loss resulted in large derivations. He was convinced that the derivations 
could be back-tracked to the pressure drop across the mud motor and the bit. These two losses can only 
be estimated by means of empirical models and are thus highly uncertain. Now he had the chance to 
record pressure loss without these two disturbing pieces of equipment installed. After lowering the 5" * 
4.127" DP (without the bit) down to the casing shoe he circulated for 45 min. to neutralize temperature 
effects. The casing, a 13 3 /s", 68 lb/ft, (ID=12.40") had its casing shoe at 4 500 mMD. The pump was a 
relatively new (volumetric efficiency = 0.96) Garden-Denver PZ-11-1600 HP triplex mud pump, 6" liner. 
One complete pump stroke delivered 15.291. At the following pump speeds, with no drill string rotation, 
he read the average stand pipe pressures (SPP) which was the average of 3 tests: 

5 spm - 13 (+/-2) bars 
20 spm - 70 (+/-3) bars 
50 spm - 100 (+/-5) bars 
100 spm - 200 (+/-10) bars 

Rheology of drilling fluid at this temperature is identified with the one in Task 2-1. Apply the Power -aw 
model. Mud density was 1.21 kg/1. The recorded pressure losses are presented in Figure 3-10. Estimate 
pressure losses and compare them with the recorded ones. 

3.10 Pressure loss. Effect of rotation 

In this task you need to use your imagination. The driller made now an additional test: At the lowest 
and the highest pump speeds [5 and 100 SPM from Task 3.9 above] he rotated the drill string at 100 
RPM for a short time and saw that the average reading changed to 8 and 248 bars respectively at the 
two selected pump speeds. Determine the effect of rotation on the mud s rheology for the given flow rate 
(the rheology is obviously dictating the annular pressure loss). Assume that the rotational movement is 
additive to the axial flow with respect to shearing effect on the fluid. When determining the shear effect, 
simply use the average rotational velocity across the annular gap. Assume also that the drilling fluid is 
the same as in task 3.9, a Power-law model with n = 0.5, K = 1.63. 
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Change or remove; too many assumptions 


200 


50 
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0 



* SPM 


0 


25 


50 


75 
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Figure 3-10: Recorded SPP (thick line), with rotation (thin line). 


3.11 Pressure loss. Nozzles. OFU 

After having drilled a 17.5" hole the 13 3 /s" casing was set and cemented at 2 300 mMD. Finally a 1214" 
hole was drilled down to the reservoir at 2 500 m. 

The drill string consisted of a drill pipe (54.276") and 100 m of drill collars (6.25 • 3.1"). The circulating 
rate was 800 GPM during drilling and the mud density was 12.9 PPG. A Fann-VG viscometer gave the 
following readings: 


Assume the mud rheology is best described through the Bingham model. OFU-equations are copied 
from Applied Drilling Engineering SPE-text book Table 4.6: 


N Re =928pvd h /ju eff 


v=q/(2A48d 2 h ) 


AP lam / Al = 


150O/ 2 22 5d 


V P i v , 



The OFU units are: p PPG 


q GPM 

A in 2 


P P S1 

C. 0.95 

d 
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a) The pressure drop through the annulus above the BHA was equal to 100 psi, but must be 
calculated along the BHA. What is the equivalent circulating density at a depth of 2 500 m? 

b) The bit had 5 nozzles, each of a diameter og 14/32 inches. Pressure loss through the surface 
pipes was 200 psi. What pressure is required from the pump when drilling at a depth of 
2 500 m?. Compare the results with bit pressure loss estimated in Si-units 

A Pbu =i 

c) Determine the pressure loss in the complete circulation system 

3.12 Swab pressure. Cling factor 

Due to the no-slip conditions on all surfaces, the mud will also cling to the drill string. The cling factor 
is used during estimation of surge & swab pressure. How would you, in a stepwise fashion, go about to 
define and estimate the cling factor? 


3.13 Swab pressure model 


Assume you are tripping out while simultaneously pumping. Your task is to start the process of derivation, 
which later, will lead to an expression of surge pressure during laminar flow. When making a drawing 


of the process, use parameters like v p (pipe), q p (pump), R w (wellbore), R p (pipe), R 0 (the point where the 
flow velocity is zero), etc, as required for your explanation. 
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4 Hydraulic program 

The exercises which are related to the hydraulic program distinguish between two different approaches 
of preparing the hydraulic program: 

a) The standard method: Each liner represents the pumps capability The hydraulic program is 
planned for one well section at a time. 

b) The extended method: All the piston sizes are treated as part of one process, and they are 
divided into two operating ranges. The hydraulic program is planned for all sections in one 
common operation. 

Both the methods are based on maximizing the ROP, which (in this book) is expressed through the 
following equation: 

ROP = A*(q/d no-i| ,)-> 

4.1 Mud pump issues 

a) Characterize a mud pump as detailed as possible with respect to 

• Effect 

• Efficiency 

b) Why are several mud pumps sometimes arranged in parallel or in series? 

c) Compare centrifugal with piston pumps 

d) Explain the term hydraulic knocking in pumps 

e) A tri-cone bit has 3 nozzles; each nozzle is 15/32 nd inch in diameter; p mud =1.3 kg/1; drilling at 
a depth of 2 500 mMD. At two pump rates (which are close to the actual operating flow rates), 
the following pressures (stand-pipe pressures) were recorded: 


q (Ipm) 

^pump v ~ ' 

P p (bar) 

2 000 

230 

670 

33 


Determine the value of and m in the expression of the parasitic pressure in an oil well. 
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4.2 Optimal nozzles? Section wise 

A tricone bit is equipped with 3-15/32" nozzles. The following data are given: 


Ppi 

P P 2 

Pbi, 

P 

r pumpe.max 

q 

Tpumpe,max 

q 

J-mm.ann 

fl-max-ann 

D 

P mud 
flopt 


= 186 bar at 2 000 1/min (p is measured at the standpipe) 

= 120 bar at 1 750 1/min 
= 1.11-1/2 pv 2 . 

= 270 bar for piston in use 
= 2 700 1/min for piston in use 

= 1 700 1/min (below this value cuttings will accumulate) 

= 2 600 1/min (above this value the wellbore adjacent the BHA will start to 
erode) 

= 3 000 m 
= 1 400 kg/m 3 

pump. max 


- 1/m 


K l D(m + 2) 


a) Determine the two constants in the parasitic pressure loss equation. 

b) Determine optimal flow rate at 3 000 m depth. 

c) What are the optimal nozzle size at 3 000 m? (in terms of x/32 inch). 


4.3 Liner selection. Section wise 

The operating data of a National 12-P-160 (this number indicates al600 HP pump) triplex pump are 
presented in this books Chapter 6 - Supportive Information. 


D = 2 500 m 
p = 1 200 kg/m 3 
K x =1.6 TO 6 
m = 1.7 

q r - 0.018 m 3 /s (d bit = 12 ^ 
q am (turb)= 0.04 m 3 /s 

a) Derive an expression of q opt and determine numerically the optimal liner at this depth? 

b) Select the most optimum liner at this depth 

c) Determine the optimal bits pressure when the 6" liner is used 

d) At what depth would you change from 6" to 5%" liners? 
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4.4 Hydraulic program. Section wise (i.e. all liners are treated as if in range I) 

a) Assume the rate of penetration (ROP) is a function of bottom hole cleaning: 


ROP = A • 


f \ 

3 _ 

\ d eJ 


«8 


Show that ROP will decrease with depth when the pump pressure is expressed through the 
equations below: 


Ppump ~ Vloss TPbit 

Pioss =K x Dq m 

Pbu = l - n \f » 2 
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b) The 12V4" section starts at the 13 3 /s" casing shoe at 1950 m MD, and is planned to reach a depth 
of 4000 mMD before setting the next casing. The following is known: 



= 

2 • 10 6 

m 

= 

1.65 

Tr, vertical 1 

= 

0.025 m 3 /s 

4r, horiz. 

= 

0.040 m 3 /s 

Tnax, vertical 

= 

0.035 m 3 /s 

Tnax, horiz. 

= 

0.045 m 3 /s 


What flow rate and liners size would you recommend through this depth interval when 
drilling either a vertical well or a horizontal well? Use the 1600 HP pump as defined in 
Supportive Information. 

4.5 Optimal parameters for BHHP. OFU. Section wise 

This exercise includes Oil Field Units, just to indicate for you how much simpler it is to work with SI units. 

The bit has 3 • 12/32" nozzles, and p mud =10 PPG while drilling at 8 200 ffMD. The following pump rates 
and pump pressures (which are within the operating flow rates q r = 240 GPM) were recorded while the 
bit was close to the bottom of the well: 


q (GPM) 

“pump v ' 

P p (Psi) 

500 

3000 

250 

800 


The pump is characterized through: 

p = 3620 psi, E = 1000 Hp 

r max r p,max r 

Determine optimal pump rate and nozzle size when applying Bit Hydraulic HP (BHHP) as the 
optimization criteria. Assume that the rate of penetration is linearly related to it. The pump volume 
efficiency is 0.9. 

BHHP = Ap bit ,( l (H) 

1714 P 

Pressure drop in OFU are ( p ( ppg ), q ( GPM),A(in 2 ), c d = 0 . 952 ) 

A Pbit = 8.311-10' 5 -p [q/(C d .A 0Z J] 2 
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4.6 Liner selection. Complete well 

The characteristics of a 1600 HP piston pump are found in Supportive Information, while the optimal 
operational functions are presented graphically in Figure 4-6. The optimal pump pressure values are in 
general the maximum ones, given for each liner size. Maximum values are the recommended values in 
the liner table, which in fact are around 85% of the absolute maximum. In Figure 4-6 the flow rates are in 
correct scale, while Revalues are only qualitative. Hydraulic parameters, the bit program and minimum 
and maximum flow rates are presented below and in Table 4-6: 


K 

K 


1 in range II 
1 in range I 


1.80 • 10 6 , m n = 1.6 
2.20 • 10 6 , m, = 1.5 


Bit diameter 

Start depth 

Permissible annular flow rate (m3/s) 

In 

m 

min 

max 

36 

0 

0.035 

0.050 

26 

100 

0.030 

0.045 

17 Vi 

500 

0.027 

0.035 

12 y 4 

1500 

0.015 

0.030 

8y 2 

3000 

0.010 

0.020 


Table 4-6: Bit program and flow rate ranges 


a) Operating range I is defined as the pump operating range of the smallest liner, range II is defined 
by the remaining liners. Derive optimal flow rate, q , in pump area II by maximizing the ROP: 


Voptn 


< £ ,), 


P -'max 


K x D(m + 2) 


1 

m +1 


b) Determine at which depth Range II stops (while drilling downwards) and at what depth the 
maximum flow rate turns into the theoretical optimum flow rate of Range I. 

c) Draw also into the graph the optimal hydraulic program (the graph is a principal drawing of 
ROP vs. q and thus not a quantitatively correct drawing). 

d) What pump rate is optimum at 2 000 m. 
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Figure 4-6: Hydraulic Data. 


4.7 Liner selection. Complete well 

A 1660 HP pump is used (see supportive info). 


a) Find at what depths the transition from operating area II to I occur. A 1660 HP pump is used 
(see supportive info. K x and m are the same for both ranges.) 

K x = 1.72 • 10 6 
m = 1.51 

b) Make a flow chart of a computer program of how to determine when to change from working 
area II to I during drilling. 
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5 Well challenges 

5.1 Filtration control 

a) How can we plan the mud composition in WBM to minimize the fluid loss through the 
filter cake? 

b) If the fluid loss shows an increasing tendency during drilling, how is it detected and how is 
the problem treated? 

c) Two sand formations of nearly equal pore pressure are encountered. Will filtrate invasion be 
greater in the sand of high permeability and high porosity, than in one with porosity? The final 
filter cake permeability in both sands is assumed to end up at around 10 3 mD. 

d) A reduction of water flow into shale is beneficial because this will reduce unwanted reaction 
between the drilling fluid s water phase and shale further away from the wall, where the water 
activity of the pore water may be different. How can the water-flow into clay be controlled? 
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5.2 Filtration control 

We want to obtain a physical picture of how far the filtrate and the particles penetrate a porous formation 
and gradually stops due to a tight filter cake. A 15" hole with open hole length of 4 000 ft is being drilled. 
The bottom 10% of the borehole length is of porous formation (defining the filter area) with a porosity 
of 15%. Assume that this porosity corresponds to the porosity of the filter paper. The filter area A of the 
filter press is 45 cm 2 (r = 3.9 cm.). 

A laboratory test of the mud showed an API water loss of 25 ml/30 min. The cumulative loss is proportional 
to the square root of time. Assume therefore that the accumulative fluid loss V is expresses as: 

V f = A • . yj = A • C • Vt 

I 

We assume the parameters defined by C are constants. 

a) To simplify the fluid loss estimation imagine that the time of drilling the well is negligible. 
Construct a plot of filtration loss vs. time. Estimate the fluid loss after 24 hours. 

b) Calculate the radius of the invaded zone after 24 hours, assuming 100% displacement of the 
pore fluid. 

c) We want to minimize the fluid loss to porous formation during overbalanced drilling, both 
with OBM and WBM. 

How would you specify the drilling fluid (focus on the part related to filtrate loss)? 
How would you follow up filtration control during the drilling phase? 

Why is filtration control important? 

The lab filtration showed an increasing tendency during drilling. Why? 

5.3 Cuttings concentration 

A horizontal section has been drilled at more or less constant ROP and flow rate. The cuttings 
concentration generated at the bit during drilling is c l = 0.02. Discuss what could be the concentration 
at these positions: 

a) At the end of the horizontal section ( = c 2 ) 

b) At the surface, when the mud is entering the return flow line (= c 3 ) 

c) What determines the cuttings bed height in the horizontal section? 

d) Why is cuttings accumulation in wellbore expansions (washouts) a problem during tripping? 

e) Mention 5 downhole problems related to poor solids control during drilling, and explain why 
or how poor solids control is the cause behind the problems. 
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5.4 Cuttings concentration 


a) What forces and mechanisms are involved when cuttings are transported in horizontal 
wellbores? 

b) Slip velocity of perfect spheres is 


F 


d cutting.g{pcutting Pmud) 


settling 


6nfd 


■m 


What is the meaning of f(c) in the given equation? Present a graphical representation of f(c) 
vs. particle concentration. 

c) Figure out with level and mentioning the forces are involved in cutting transportation in high 
deviation wellbore. 

d) Cleaning of horizontal wells is a challenge. Your task is to: 

• Explain the principles of how cleaning works and which processes and parameters 
are involved. 

• Why does the drill string RPM needs to be > 120 RPM before the cleaning process 
become really efficient? 
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5.5 Density control 

a) Define a weighted mud system (as opposed to an un-weighted)? Explain how to clean 
weighted muds. 

b) Derive a simple formula of necessary volume increase, AV dd , involving weight material with 
density p a dd (4.3 kg/1) to increase the density from pi to P 2 . Use this formula to estimate how 
much mass of barite must be added to increase mud density from 1.3 to 1.4 kg/1. Original mud 
volume was 60 m 3 . 

c) 1 m 3 of mud has a density of 1.5 kg/1. Adjust the mud density to 1.72 by adding 100 1 mud of 
density 1.8, 40 kg Bentonite (to adjust rheology) of density 2.3 kg/1 and barite of density 4.2 
kg/1. Find how much Barite of density 4.3 kg/1 is needed to obtain a mud density of 1.72 kg/1, 
all ingrediences added together simultaneously. 

d) Different water based muds defined below are stored in three different mud pits. All 3 pits 
should be mixed into one tank and water added until the density becomes 1.55. How much 
volume of water must be added? 

Vj = 10 m 3 , p i = 1.5 kg/1 
V 2 = 20 m 3 , p 2 = 1.6 kg/1 
V 3 = 3 m 3 , p 3 = 1.9 kg/1 

5.6 Density control 

Sometimes the mud viscosity increases unintentionally due to accumulation of fines in the mud. These 
fines are referred to as Low Gravity Solids Content (LGSC). The fines are too fine to be removed by the 
cleaning equipment. Typical density of LGS is 2.4 kg/1. They are inert, but builds viscosity since particle 
size is small (< 5 p). Their unwanted effect can be reduced by diluting the mud with water. Here follows 
3 examples: 

a) The mud volume is 100 m 3 with a density 1.8 kg/1. The fraction of low gravity solids is too 
high, 5 weight %, and has to be decreased to 3% by water addition. Calculate the mud volume 
to be discarded and the amounts of fresh water and Barite that should be added. The original 
volume and density has to be unchanged. 

b) A tank containing 90 m 3 of mud has a density of 1.6 kg/1 and should be increased to 1.7 kg/1. 
The volume fraction of low-gravity solids must first be reduced from 0.055 to 0.030 by water 
dilution. It is required that you first discard a part of the original mud volume, so that after 
adding of water the volume of the mud is 90 m 3 before the barite is added. How much barite 
must be added, and what exactly is the new LGSC? 

c) Calculate the volumes of old mud (< 10 m 3 ) and barite that has to be mixed in order to fill a 
10 m 3 large pit with mud which must balance a pore pressure of 410 bar in a depth of 3000 m. 
Barite has a density of 4.3 kg/1. The density of the old mud is 1.2 kg/1. 
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5.7 ECD. Barite 

Explain the reasons behind and suggest potential solution to the following problem: 

This well was drilled with WBM, weighted by Barite to 15 PPG. While POOH to change the \2 X A" bit, 
the driller experienced no problems. When GIH with the new bit some weight reductions (took weight) 
were experienced in the build-up zone. It took around 5 h from the bit left the bottom of the well till it 
returned. After the bit reached the bottom, the well was circulated for some time, and the returning mud 
behaved strangely. The mud weight, which originally was 15 PPG displayed an initial sharp decrease, 
then increased again as shown in Figure 5-3, before finally stabilizing at 15 PPG. 



5.8 ECD. Fluid and flow 

A 17Vi" hole was drilled from the 20" casing shoe at 1100 mTVD to 2 100 mTVD. The bottom hole 
assembly consisted of 120 m of 9 V 2 " Drill Collars (DC). A 5 V 2 " drillpipe (DP) was used. 

Capacities: 


I 7 V 2 " open hole capacity: 

155.2 1/m 

DC / Open hole capacity: 

109.4 1/m 

DP Open hole capacity: 

139.2 1/m 

DP / Casing capacity: 

161.8 1/m 

9 V 2 " DC / capacity: 

4.56 1/m 

5 V 2 " DP capacity 

10.77 1/m 


Mud Parameters: 


• Mud density: 


1.25 kg/1 

• Rheology: 

600 / 300 rpm: 

51.7 / 30.6 Pa 


200 / 100 rpm: 

22 / 12 Pa 


6/3 rpm: 

3/4 Pa 


Gel: 10s / 10 min: 

5 / 13 Pa 
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a) Prior to drilling, the hole was circulated at 3 500 1/min. What is the annular pressure loss, and 
what is the corresponding ECD? 

b) The drilling commenced from 2 100 m. The same rheology and flow rate was applied. At 2 
300 m the average drilling rate was 50 m/hr. Formation bulk density was 2.4 kg/liter. What 
is the ECD in this situation? Transport ratio is 0.75. At the casing shoe at 1 800 m TVD the 
formation fracture pressure was 228 bar. Check if everything is OK. 

c) A new mud was being prepared for the 8Vi" section. The 13 3 / 8 " csg shoe was located at 15 000 
ft vertical depth. While drilling at 16 500 ft the well started losing mud and it was decided 
to lower the MW from 15 to 14 PPG. The well had very narrow pressure window, and the 
equivalent pore pressure gradient at this depth was 13.5 PPG. The mud was mixed to 14 PPG 
with an effective viscosity of 40 and 30 cP at 600 and 300 RPM respectively, at an average surface 
mud temperature of 40°C. After the new mud was circulated, the pump was shut off, and a 
flow-check indicated that the well was dead (no influx). While repairing the power swivel the 
well started to flow by itself, and soon afterwards the kicking well had to be shut in to avoid 
a complete unloading. 

How would you go about to estimate the pressure profile of an initially cold, static fluid column 
as a function of time. You are asked to present the governing equations of heat transfer in a 
well on differential form, and work out a flow sheet of how to numerically solve this task. 
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5.9 Water activity 

a) Define water activity, A w and how to determine A w in a a) salt water solutions and in b) pore 
water in shale? 

b) Explain why water activity is a function of water salinity? 

c) Explain why high water activity causes clay swelling problems? 

d) How do you prevent clay swelling problem? 

e) How does water activity of the water phase in OBM influence wellbore stability? 

5.10 Shale stability 

In order to avoid that water enters and causes the shale to swell, the activity of the water in the mud 
(including the water phase in oil based mud) and in the shale must be equal. 

The activity of the water in the water phase in shale cuttings is measured in the field using an electro 
hygrometer. The probe of the electro hygrometer is placed in the vapor above the sample being tested. 
The electrical resistance of the probe is sensitive to the amount of water vapor present. Since the test 
always is conducted at atmospheric pressure, the water vapor pressure is directly proportional to the 
volume fraction of water in the air/water vapor mixture. The instrument is normally calibrated with 
saturated solutions of known activity shown in Table 5-6. 


Salt 

Activity 

ZnCI 2 

0.10 

CaCI 2 

0.30 

MgCI 2 

0.33 

Ca(N0 3 ) 2 

0.51 

NaCl 

0.75 

(nh 4 ) 2 so 4 

0.80 

Pure water 

1.00 


Table 5-6: Saturated solutions of different salts and its vapor's water activity 

Sodium chloride and calcium chloride are the salts mostly used to alter the activity of the water in the 
mud. Calcium chloride is quite soluble, allowing the activity to be varied over a wide range. In addition, 
it is relatively inexpensive. The resulting water activity for various concentrations of NaCl and CaCl 2 are 
shown in Fig. 5-6. 

a) The activity of a sample of shale cuttings drilled with OBM (no foreign fluid invasion) is 
determined to be 0.69 by an electro hygrometer. Determine the concentration of calcium 
chloride needed in the water phase of the mud in order to have the activity of the mud equal 
to the activity of the shale. 


39 


Download free eBooks at bookboon.com 




Exercises in Drilling Fluid Engineering 
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Figure 5-6: Water activity in calcium chloride and sodium chloride at room temperature. 


b) A core is taken from a swelling formation. Can you retrieve any useful information from its 
specific weight, useful with respect to avoid swelling while drilling through it? 

c) Explain why wellbores and cuttings stability is so much better when applying OBM instead of 
WBM. As part of the answer, please explain the principal function of the two different surface 
active additives that are always added to Oil based mud. 

5.11 Shale stability 

a) When drilling into swelling clay, problems like sloughing (soft) shale and stuck pipe may occur. 
Explain how/why this can be avoided by means of the proper oil based drilling fluid and specify 
the ingredients in the drilling fluid. 

b) Why is KOH preferred over NaOH? 

c) What significance does the K + concentration have for the shale? 

d) Can wate flow through shale be controlled? 

e) Explain the principal function of the two different surface active additives that are always 
added to Oil based mud. 

f) How is the salt concentration in the water phase, which is added to Oil Based Mud, determined? 


While drilling in the 8 V 2 " section, at a depth of 1 500 mTVD / 6 000 mMD, in overbalance, the ECD 
will fluctuate and at times be high in this long well. Previous experience from that area indicates that 
instable, swellable shale will be penetrated. Your task now is the following: 

g) Define what wellbore stability-related processes may take place in the shale while drilling 
through it with WBM. 

h) Which type of inhibitive mud will you suggest in order to maximize wellbore stability? Explain 
how this mud type will affect the wellbore. 

i) Does fluctuating ECD have any implications for the stability of the wellbore? 
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5.12 Wellbore problem 

a) What are the dominating mechanisms or factors leading to mechanically stuck pipe. Explain 
the mechanisms of differential sticking in porous/permeable formation. 

b) What are the consequences of stuck and how do you suggest combating the problem? 

c) What is the most likely stuck pipe mechanism while drilling in salt formations? In the case 
of presence of halite type salts, what type of mud would you select for safe drilling in such 
salt section. 

d) Wellbore breathing (ballooning) and Seepage losses. Include the headings; Definition, 
Explanation, Detection; Repair activity. Discuss the two phenomena. 
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6 Supportive Information 


6.1 Pump (National 12-P-160) and hydraulic program data 


Line size 

in 

5V2 

5% 

6 

6 14 

6 Vi 

6% 

7 

7 V 2 

Discharge 

Psi 

5555 

5085 

4670 

4305 

3980 

3690 

3430 

3200 

pressure 

10 5 -Pa 

383.0 

350.6 

322.0 

296.8 

274.4 

254.4 

236.5 

220.2 

Pump rated at 

GPM 

444 

486 

529 

574 

621 

669 

720 

772 

120 spm 

m 3 /s 

0.0280 

0.0307 

.0334 

.0362 

.0392 

.0422 

.0454 

.0482 

HP 


1439.0 

1441.8 

1441.3 

1441.7 

1442.0 

1440.3 

1440.8 

1441 


Power of Efficiency = 1441.7 • 745.7 = 1.0748 • 10 6 or 322 • 10 5 • 0.0334 = 1.0755 • 10 6 (watt) 


ROP,u. =C\q!d,f 

ROP mHr =C\BHHPf 
P p = q-P p (Power) 


P P = A Pbit + A Pd 
A Pd = K i D • °r 
p bit =l.U^pv 2 


l optj-q / d n 


2 P t 


(m+2)K x D 


q, 


opt u -q/d n ' 


( £ ,)„ 


K l D(m+ 2) 


1 

m+\ 


BHHP = Ap bit ■ q 
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6.2 Pressure loss equations 


Newtonian fluid 


Bingham model 


Power law model 


Lam/pipe 


Lam/annulus 


Turb/pipe/ann 


Ap, 

A Pa 


32 v juL 
p d 2 
48 v p L 

{d 0 -d t f 


32 p-L-v 16 Lt 

Ap n = --+-- 

p d 2 3d 

48 p p L-v 6 Lt 
A p a = -— +-~ 


, 8v 3n + 1 \ L 

Ap=4K\ — •- • - 

p 1 d 4 n d 


Ap = 


d\ 2 


Ap 


( d 0 ~d i ) z d o -d i 
0.073 • v 1 ' 8 • p° p f • L 


A Pa = 4 K 


12v 2n +1 


d„ - d. 3 n d„ - d. 


4 L 1 


Ap = a ■ N Re - — - - pv 


d h r 


a = (log n + 3.93)/50 
b = (1.75- log n)/l 


Eff. vise, pipe 
Eff. vise, ann 

Shear-r. pipe 

Shear-r. ann. 


Meff=T/y 


Veff= Vy 

. sv 

Y = v 


Veff = Vpl + 

^eff = + 


6v 

T o ( d o ~ d , ) 




-( 


8v 3/7+1 


8v 


Metr =\ 


d 4/2 
12v 2/z+l \ Kd 


Kd 

8v 


3/2 


12v 


7 = 


12F 
d . - d. 


. 8v r 
y= — + — 
d 3 u p/ 

YLv r 

r =-—- + - 


In +V 


\d 4 n j 

\2v 2n +1 


d-d. 2a 

n i * 


Pi 


x d -d. 3n 

\ O l 


General 
N • 

^ 1 re,pipe 

General N rP 


N Re = 




cT -v 2 ~” -p 

V( s'" 1 ) 

dp v 2 ~"-p 

Prl 12- 1 ) 


x,-*- 




3 / 2+1 

4/2 


2 / 2+1 

3/2 


Fanning f lam = 16/N re 


Fanning f lam = 24/N re 
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Continuity equation: _ ^R_ - y. (py) 

dt } 

Microscopic Cylindrical coordinates - — = — — ( v prv r ) + — — (pv 0 ) + — (pw) 

dt r 8r r r SO 8z 

d\vpdV 

Macroscopic —-= -A pvA = p x v x A x - p 2 v 2 A 2 

dt 

Dv 

Momentum equation p — = pg - Vp - V • r 


Microscopic Cylindrical coordinates (only the r-component) 


dvr dv r Vn dv „ dv, 

—- + v —- + — —- + v —- 
dt dr r d 6 dz 


dp 

= —— + ju 
dz 


dr 


1 d / x 
r or y 


+ - 


1 0 2 v. d 2 v. 


r 5^ dz^ 


Macroscopic — f pdV = p x v x A 1 - p 2 v 2 A 2 + p x A x - p 2 A 2 - F + Mg- 
dt J 

The steady state, one dimensional pipe flow form is: p x A x — p 2 A 2 —F — Mg sin# 
Energy equation 
Macroscopic 


y -2 

p v 

— + — + z 


r 2 g 


+ A 


pump 


f -2 

£> V 

— +-+ Z 


y 


/ 2 g 


+ h 


friction 


y out 


6.3 


Conversion factors and formulas: 


1 Hp: 

Shear stress: 
Shear stress: 


l si ■ 


745.7 W 

x 0FU = 0 -1.06 (Fann VG readings = 0) 
1 lb/ 100ft 2 (OFU) = 0.4788 Pa (SI) 
t ofu -0.4788) 


Shear rate: 

y(s~ x )=RPM- 1.703 

1 inch: 

0.0254 m 

1 bar: 

10 5 Pa 

lcP: 

10 3 Pas 

P (effect): 

q • p (Watt) 

^"sphere* 

4 . Jtr 2 

V h : 

sphere 

4/3 • Jtr 3 
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i 


Content: 

1 . 

2 . 

3. 

4. 

5. 


1.5 


a) 


b) 


c) 


d) 


Solutions to exercises in drilling 
fluid engineering 


Fluid Properties 
Rheological models 
Drilling fluid dynamics 
Bit hydraulics 
Wellbore challenges 

Rheology control 

The mean size of Barite is typically 20 mm. They will be uniformly dispersed in the drilling 
fluid and they will lead to increased viscosity. PV and p eff will increase. YP will be unaffected 
or decrease, since the smallest Barite particles will behave as physical dispersants. 

In order to: 

Suppress Ca ++ from dissolving. 

Keep anionic colloidal particles dispersed. 

Suppress corrosion, H 2 S- and C0 2 - attack. 

Polymers are polymerized monomers. They can be of organic origin or be manufactured 
synthetically. Polymers have high molecular weight and come mostly as charged particles -> 
they will bind water molecules -> increase the hydrodynamic volume -> influence viscosity 
and filter behaviour. Some types of polymers can attract or bind charged particles (clay), and 
therefore influence the swelling process and contribute to selective flocculation (clear water 
drilling). The purpose of polymers are mainly to: 

• increase fluid viscosity 

• reduce fluid viscosity in turbulent flow (drag reducer) 

• control flocculation 

• improve both the filter cake itself and the fluid loss through the filter cake 
The description of the three expressions are: 

• Pseudo plastic fluids exhibit altered apparent viscosity whenever the shear rate is changed. 

• The fluid displays a reduction in viscosity over time at constant shear rate as indicated 
in Figure 1-5. 

• Rheopectic fluids are rare. They exhibit increased shear stress at increasing shear time 
(at constant shear rate). 
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Figure 1-5: Behavior of thixotropic fluids after stillstand (left) and at constant shear rates. 


e) Clay particles have a static, negative surface charge localized on the particle edges, but with 
weakly positively sites on the surface of the platelets. This triggers water! When studying the 
repulsive - attractive forces, it is experienced that at low salt concentration or high colloidal 
concentration, a very slow flocculation will take place. When edges come sufficiently close to 
the surfaces of other Montmorilonite particles, they join. 

1.2 Rheology control 

a) Bentonite particles attract water molecules (dipoles) in hundreds of layers onto each of its 
charged surface. Swelled Na-Bentonite sheets separate readily when exposed to shear forces. 
Salt will reduce the charges on the Bentonite surfaces. In addition the salt (ions) will bind much 
of the water (reduce the water activity). 
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b) At still-stand many layers of water molecules are attached to the Bentonite surface, and particles 
flocculate slowly -> viscosity increases. When the dispersion is stirred or pumped, the resulting 
shear stress will remove attached water layers and break up the flocks. Dispersed Bentonite 
with many free water molecules, torn off due to high shear stress, has lower viscosity than at 
still stand. 

c) Drilling in sediments of high yield (premium) clay produces the following clay concentration: 
Fraction of solids by weight =/ 


Qclay p c i a y _ ROP • Afrit' p c i a y 

Qpump'Pmud Ppump'Pmud 


0.0097-0.34-2400 
0.05 -1 025 


0.154 


ROP = 35 m/t = 35 / (60 • 60) = 0.0097 m/s 
qpump = 3 000 1/min = 3 000 / (1 000 • 60) = 0.05 m 3 /s 
Abit = tt/ 4- (dbit 2 ) = ti/ 4 • 0.662 2 - 0.34 m 2 
dbit = 26” = 26 • 0.0254 = 0.662 m 


Drilled clay contributes to a weight increase of 15.4% and produces a viscosity of 42 cP, which 
is above the required viscosity of 15cP. 


d) The yield of Bentonite is how many m 3 of mud of 15 cP which the amount of one ton Bentonite 
is able to produce. From Figure 1-6 one ton of Bentonite will produce 6 weigth% = 16.666 m 3 
of 15 cP mud. 

e) From Figure 1-6 we see that a viscosity of max. 50 cP will correspond to 7.5 w% of Bentonite. 

P b en,oni,e= 2 ^kg/m i 


To find the density of the mixture of 100 kg mud =^> 7.5 kg Bentonite + 92.5 kg water, we need 
to find the volume: 


Volume = V, 


bentonite' 


V = 

water 


1.5 kg 


_ m _ 100 kg 

Pmud ~ v ~ 0.0955 m 3 


_ + 92,5 kg 

2400 kg/m 3 1000 kg/m : 

— 1 047 kg/m 3 


= 0.0955 m 5 


This shows that mud density cannot be increased much higher than 1 047 kg/m 3 by the addition 
of Bentonite. 


Flocculation 

a) When two Bentonite flakes are sufficiently near each other they are electrostatically attracted 
and will join edge to surface. 

b) The cement is not by far hardened and contains a high concentration of lime (CaOH). Lime 
dissociates and one Ca ++ ion can crosslink two charged clay platelets. This binding cannot be 
broken by hydraulic means (shear stress), and the flocculated mud has to be dumped after 
having been circulated to the surface. 
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c) The factors which enhance flocculation: 

• High concentration of charged colloids (Bentonite or anionic polymer) 

• High concentration of Ca ++ ions (from dissolved salt and carbonates) 

• High Temperature (Brownian movements of the water molecules contribute to 
bringing the colloidal particles more frequently in contact with each other) 



-► 


d) WBM behaves shear thinning: 

• Amount of water dipoles attached to each charged colloidal is a function of the shear 
force. High flow rate-> highe shear rate -> water molecules are sheared off and become 
free -> lower viscosity of the suspension 

• Weak flocculation occurs, especially at low sehar. At high shear the floccs break up -> 
lower yeld point 

This behavior was presented during lectures. 

e) The shearing action explained in e) above also is time dependent. Layer by layer are sheared 
off in concedutive order (presented in lectures). 

f) The surface of caly colloids are positively charged by the loosely bonded Na-ions and represents 
therefore a weak bond compared to the electrostatic bond of double valence cations, mainly 
Ca ++ , to the strolngly negative charged colloids at their edges. 

Mud contamination 

a) Contaminations are dissolved salt and chalk; drilled through cement etc. all produce Ca ++ 
ions which lead to flocculation. Also fines are pollutant. At the surface the mud can be treated 
with thinners. Low gravity (2.4 kg/1) solids content (LGSC) is reduced by running the mud 
through centrifuges. 

b) Rheology monitoring. First we plot the evolution of the rheology vs. time, as shown in 
Figure l-8a. 
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Then to evaluate the mentioned changes of the Bingham field rheology-model is perfect. We estimate 
the Bingham constants: 


= 42 

- 28 

= 14 cP 

= 41 

- 27 

= 14 cP 

= 68 

- 54 

= 14 cP 

= 31 

- 17 

= 14 cP 


We observe that the mud is dispersed at first. Salt will dissolve and Ca ++ contaminates the mud. We now 
see that the viscosity is constant, while the yield point (YP) has changed like this: 


YP 

= 28-14 

= 14 


= 27-14 

= 13 


= 54-14 

= 40 


= 17-14 

= 3 
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The PV is constant (constant inner friction means solids content is constant). A YP is caused by higher 
attractive forces, caused by drilling into a salty formation. Initially, low concentration of Ca ++ ions leads 
to strong cross binding (flocculation). Later, when the Ca ++ concentration has been high for a while the 
Ca ++ ions have exchanged the Nations in the Bentonite plates, leading to aggregation of clay platelets. 
A high concentration of Ca ions will over time lead to lower shear stress than originally due to cation 
exchange. The development is seen in Figure l-8b. 



Countermeasures: Add dispersants or use Gyp mud or use a high pH level in the mud when 
contaminations are expected. This is how to avoid flocculation and aggregation. 

c) The edge-to-face flocculation is weak electrostatic forces and is easily broken when sheared. 
The Ca ++ ion binds two clay platelets in a much stronger grip, and is not easily broken as 
indicated in Figure l-8c. 



Figure 1-8c: Rheological response at the surface when drilling into contaminants. 


1.5 Flocculation 

a) Consider the force balance between shear and gravity acting on a particle: 

A p ■ t y iTip ' g Vp (pp pfi ) g 

1 

4 n dp ■ T y = —n dp (p p - p ft ) g 
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Solve for yield point 

cL , . 0.005 „ 

T y = Y\ ( pp ~ pfl 1 9 = 24 (2 300 “ 1100) 981 = 2 ' 45 Pa 

Now solve the same balance of force while particles are settling, with respect to particle diameter: 

d p = x y • 24 / [(p p - p fl ) g] = 15 Pa * 24 / [(2300 - 1100) -g] = 0.0304 = 30.4 mm 

b) The gel structure inside the drill string could have broken either by just pumping, by hitting the 
pipe with a hammer, or by inserting a heavy pill in the upper part of the drill string (U-tube 
the level in the drill string downwards). Consider the force balance between shear force and 
gravity acting on the fluid column inside the pipe: 

Afluid' 'Ey — pg ' Ah 

7rdpi pe ■ L ■ T y = pg ■ Ah 

Solve for Ah, the additional hoisting height before gravity breaks the gel 

Ah = n ■ dpipe- L ■ ^ = n ■ 4.127 ■ 0.0254 ■ 2100 ■ 14.4 /(1100 ■ 9.81) = 0.3 m 


1.6 Fluid additives 

a) 


Material 

Definition 

Relevance 

Anhydrite 

CaS0 4 = sedimentary salt 

dissolves in water -> Ca ++ 

caustic soda 

NaOH, Sodium hydroxide = lye 

adjust pH (dispersant) 

Ca S0 4 

anhydrite 

contaminator, leads to flocculation 

Gypsum 

Ca S0 4 • H 2 0 

Sedimentary salt 

Restructure from Na- to Ca-Bentonite; shale control 

CEC 

Cation Exchange Capacity 

ability of reactive clay to exchange cations 

CMC 

natural polymer 

filter loss/viscosifyers 

Chalk 

CaC0 3 = sedimentary salt 

dissolves in water -> Ca ++ 

Colloid 

particles of size <2\i, clay, silt 

build viscosity 

Deflocculator 

dispersion of colloids 

reduce/neutralize electrostatic attraction 

Dispergator 

spreading of colloids 

thinner 

HEC 

natural polymer 

viscosifyer, filtrate reducer, deflocculant. 

Lignite 

anionic polymer 

Thinner, filtrate reducer 

Lignosulphonate 

anionic polymer 

thinner (dispergator) 

Dispergator 

thinner 

hinder particles to coalesce 

MBT 

Methylene Blue Test 

determines amount of reactive clay in mud. Clays 
present in the mud will adsorb MB 
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PAC 

natural polymer 

filter reducer, shale control 

PH PA 

synthetic polymer 

shale control, Bentonite extender 

Prehydrated 

Bentonite hydrated in water 


Salt 

Ionic compound, NaCI 

dissolves in water ->Na + 

SAPP 

anionic polymer 

thinner (dispergator) 

Starch 

_"_ 


Xanthan 

natural polymer 

viscosifyer 


b) A drag reducer consists of long-chained, neutral polymers added in small amounts to the 
viscous fluid. In a turbulent flow regime, the fluid molecules move in a random manner, causing 
much of the energy applied to them to be wasted as eddy currents (and corresponding high 
pressure loss). Due to shear stress the polymers are stretched out where the shear is high (close 
to the wall). In outstretched state the diffusion is higher due to lower effective surface area 
and thus lower resistance to movements. The polymers diffuse more in the direction of the 
pipe centre since the wall hinders diffusion in the opposite direction. The concentration will 
therefore become low after some flowing distance. The small concentration of stretched out 
long polymers along the wall will tend to turn the turbulent flow into laminar or streamlined. 

c) The higher concentration, the more K + is exchanged with Na + shale, the less shale swelling is 
seen. K+ is geometrically suitable and leads to high platelet attraction. 
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d) Figure 5-10 shows both how a Ca attacks is changing the dispersed mud, but also how Gyp 
mud is made. Bentonite is composed mainly of the mineral called Montmorilonite. It is an 
aluminium silicate in which the silicon-oxygen sandwich is a layer of aluminium hydroxide. 
Some of the aluminium in the lattice structure is, during creation, replaced with ions of a lower 
charge (electron valence) such as magnesium (from 3 to 2). The substitution , without any other 
change of structure, creates surface negative charges in the lattice. 

The negative charge sites created in the clay sheet are partly balanced by close association of 
positively charged ions (cation), normally Na + . The creation of these charged ions and charged 
clay surfaces creates very strong attractive forces for polar water molecules that readily force 
themselves in between the unit layers. 

Gypsum is added to prevent shale swelling. When lime or gypsum is added to Bentonite-treated 
mud, sodium Montmorilonite will convert to calcium Montmorilonite, which first produces 
flocculation and eventually aggregation of the Montmorilonite (see Figure 5-10). Caustic soda 
is added for pH control to suppress further dissolution of Ca ++ . Lignosulphate is added for 
deflocculating. CMC may be added for fluid control. Then we get gyp mud. 

If a lot of swellable clays have to be drilled through, apply gypsum mud; Pre-hydrate Na + - 
Bentonite Add gypsum to the mud. Gypsum (Ca S0 4 ) dissolves and Ca ++ ions are formed. 
First reaction causes flocculation. Over time Ca ++ replaces Na + in Bentonite platelets, to make 
Ca ++ - Bentonite and they aggregate and viscosity goes down. But from now on, there is an 
access of Ca ++ in the mud (and in filtrate). No new flocculation will occur, and swellable clay 
does not swell much. 



Stable Na + Montmorilonite Ha* partially exchanged Face to edge and face 

dispersion by Ca** to face Flocculation 


Ion-exchange is dependent on Ca** concentration. Low concentration 4 Partly cation exchange, high concentration 
4 ccomplete cation exchange. Partly exchange leads to flocculation tendency. 



Yield point (and effective viscosity} increases when low concentration of 
salt contamination. Effective viscosity indicated by dotted line. 

Figure 5-10. Clay behavior. 
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1.7 Fluid additives 

The pH value of WBM are held at 9-11 to increase the solubility of anionic compounds. 

a) The change in OH concentration required to increase the pH from a particular value to another 
is given by: 

A [OH ] = A [OH ] „ - A [OH ] „„„ 

A [OH ] = 10 ln ' 14) - 10 ?: " " = 0.001 - 0.000000316 = 0.00099 mol/1 

b) Mass of NaOH in g/1 = concentration molecular weight (of Caustic Soda) 

= 0.00099 moles/1 (23 + 16+1) g/moles = 0.0198 g /I 

c) K + is geometrically suiting in between the Na-Montmorilonite platelets, and leads to higher 
platelet attraction (low swelling) than for Na-Montmorilonite. Recall in your chemistry or 
periodic table (Group 1) that potassium (atom # 20) is more reactive than sodium (atom 
#11). Why? The atomic radius of potassium is larger than that of sodium. Therefore, the single 
valence electron that exists for all alkali metal is located further away from the nucleus for 
potassium than sodium. This results in less energy required to remove that valence electron 
from potassium than from sodium, leading to increased reactivity. This trend continues as you 
move down Group 1 on the periodic table; i.e. Rubidium is more reactive than K. 

d) The mud must have sufficient viscosity & velocity to exceed the settling velocity of the cuttings. 

e) The mud must be thixotropic, i.e. must have the ability to gel when stationary or at low laminar 
flow, but becomes less viscous during circulation or when rigorously shaken by the shale shaker. 

f) A large proportion of the mechanical energy (in the form of WOB and rotation) and hydraulic 
energy will generate friction and is dissipated as heat. Drilling fluid, which has a high heat 
capacity, is voluminous (typically 100 m 3 ) and absorbs this heat and allows the drill bit and 
the rest of the drill string not to be intensely heated. 

g) The drilling fluid prevents blowouts by providing a hydrostatic pressure at least greater than 
the formation pressure. 
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2 Rheological models 

2.1 Bingham / Power-law 

a) Field method: n = PV = 0 m - 0 3OO = 106.4 - 75.0 = 31.4 cP 

YP = 6» 300 -PV = 75 - 31.4 = 43.6 lb/100 ft 2 = 43.6 • 1.06 • 0.4788 = 
22.1 Pa 

Standard method: u , = — = 54 0-381 = 0.0311 Pas = 31.1 cP 

p Ap 511 

YP = t 30 o — Hpi • Y = 38.1 - 0.0311 • 511 = 22.2 Pa 
liNewton = ~ = ^ = 0.053 P(IS 

7600 1022 

n = 3.32 log (54/38.1) = 0.50 
K = 38.1 / 511° 5 = 1.69 
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Exercises in Drilling Fluid Engineering 


Rheological models 


b) 100 PPM = 100 • 1.703 = 170 s 1 

At a shear rate of 170 s 1 the shear stress is seen in Table 2-1 to be 21.5 Pa. To answer the 
question please estimate the shear stress at this shear rate and compare. 


^B-H,170 

^P-L,170 


= 27.4 Pa 
= 21.5 Pa 


P-L is obviously best. 


c) To answer the question we need to know at which shear rate to compare the two models. Given 
data gives us v = q/A = [1 000 / (1 000- 60)] /- • 0.1 2 = 2.12 m/s. 

4 


This results in shear rates for the three models, taken from Chapter 6.2 in this exercise book. 
Pressure loss equations, equal to: 


2.12 


if newt — 8 ^ — 8 — 170 S 


.-1 


W=8 Z + i£= 170 + 


0 v 3n+l 

Vpj =8-- 

' d 4 n 


3 \i p i ' (3-0.0311) 

8-2.12 3-0.5+1 


= 407 s" 


o.i 


4-0.5 


= 212 s" 1 - 


Checking the three models theoretically, at the estimated shear rates: 

T newt = M 7 = 0.0601 ■ 170 = 10.2Pa 

iBingh = t 0 + /ipi ■ y ‘=22.1 + 0.0311 • 407 = 34.7 Pa 

tpl = K-y n = 1.63 ■ 212 0,5 = 24.4 Pa 


Then the recorded; From the graph we read the following shear stress 


Model 

Shear rate 

Read 

Theoretically 

A 

Newton 

170 

21.5 

10.2 

-11.3 

P-L 

212 

23.4 

24.4 

+1.0 

Bingham 

407 

33.2 

34.9 

+1.7 


From this investigation we may conclude that P-L and Bingham model predict the rheology rather well 
at this flow rate. 
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Exercises in Drilling Fluid Engineering 


Rheological models 


2.2 Bingham/Power-law 

a) The flow curve (x vs. y) data show obviously closer nearness to the Power-low model, especially 
for low shear rates. The Bingham model reads about 12 Pa higher than the measured data. To 
verify we need model constants: 


log(71.05 / 49.75) log 1.428 

n = -= —---— = 0.51 

log(1022 / 511) log 2 



71.05 

1022 0 - 51 


2.07 Pas~ n 


P v i = 06oo - 03 oo = 140 - 98 = 42cP, r 0 = 0 3 oo “ P v i = 98 - 42 = 56 

= 56 ■ lb/100ft 2 
= 56 ■ 1.06 ■ 0.478 = 28.42 Pa 


Verification at 100 RPM (170 s' 1 )* where we read shear stress to be 27.39 


Tioo ,PL = K-y n = 2.07 ■ 170 0 ' 51 = 28.41 Pa 

Tioo ,bh = T o + Ppi'Y = 35.6 Pa = 28.42 4- 0.042 • 170 = 35.56 Pa 


We observe that the Bingham over predict as expected, while P-L fits well. 

b) This is a time-effect which occurs while performing rheological readings in the lab. The 
time being spent while performing the test at low speed gives the colloidal clay particles the 
possibility to flocculate (edge to face). The flocculation process is much more pronounced at 
low shear rates. 


2.3 Bingham/Power-law. Regression 

a) All model constants are necessary for plotting of the three models: 
Newtonian model: We initially select the 300 rpm reading 


M 


T 

7 


42 • 0.4788 
300 • 1.703 


= 0.039 Pas 


Bingham model (Field procedure): 

Ppi = 06oo _ 03 oo = 60.4 - 39.6 = 20.8 cP = 0.021 Pas 

z y = T300 “ Ppi = 39 - 6 _ 20 - 8 = 18 - 8 = 18 - 8 ‘ i- 06 '0.4788 = 9.6 Pa 

Bingham standard procedure: 

[i pl = ^ = (30.6 - 20.2) / 511 = 0.021 Pas 

Ty = t 300 - p p i ■ y = 20.2 - 0.021 • 511 = 9.5 Pa 

Both procedures gave the same answers. This assumes that the factor 1.06 is taken into 
consideration in standard procedure. 
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Rheological models 


Exponent model: 

tOg^QQ- l0 g r 300 

iog/eoo- log/aoo 


log- 


n = 


logZ«00 


log 


64 

42 


log 


64 

42 


600 • 1.703 


log 


log 2 


= 0.60 


K = — =- 


64 • 0.4788 


r 3 oo 300 • 1.703 

= 0.454 Pas n 


y n (600 1.703) 06 

To plot the three curves (see Figure 2-3.1), x is calculated at 50, 100 and 800 rpm: 


Bingham: 

T 

= 

9.5 + 0.021 • 50 • 1.703 = 11.3 Pa 


T 

= 

9.5 + 0.021 • 800 • 1.703 = 37.8 Pa 

Power-law: 

T 

= 

0.454 • (50 • 1.703)° 61 = 6.8 Pa 


T 

= 

0.454 (100 • 1.703) 061 = 10.4 Pa 


T 

= 

0.455 (800 • 1.703) 061 = 36.5 Pa 
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Exercises in Drilling Fluid Engineering 


Rheological models 



Readings on flow curve: 


50 rpm = 10.0 Pa 

100 rpm (already measured) 

800 rpm = 36.5 Pa (extrapolated) 


From this we see that the Bingham model suits the mud is flow curve best at 50 s' 1 , since x 50BM is 
much closer to the real data than x 50p L , while at higher shear rates Power law is the better one. 

b) Plug flow will occur at any laminar flow since shear stress is low in the middle of the pipe. The 
Bingham model predicts an YP = 9.1 Pa. This is a theoretical value, but given time at steady-state 
laminar flow, gel strength will develop, in fact approaching the level of the YP. Shear stresses 
below this level will lead to gel (flocculation) as shown in Figure 2-3.2 



Figure 2-3.2: Plug flow of a Binghamian fluid. 


2.4 Effective viscosity 

a) The effective viscosity (Newtonian model): 


For 600 rpm: 
For 3 rpm: 


43-0.4788 -1.06 
600-1.703 


= 0.021 Pas = 21.2 cP 


M = - = - 

r 


3.2 • 0.4788 • 1.06 
3-1.703 


= 0.318 Pas = 318 cP 
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Exercises in Drilling Fluid Engineering 


Rheological models 


These and more data points result in Figure 2-4: 



Figure 2-4: Viscosity for a non-Newtonian fluid. 


b) When determining the effective viscosity we assume the fluid is closest to a Power law fluid: 


2 data points-procedure 

log t 2 /t 1 _ log 43/30 
log 72 M log 2 


n = - 


„ t 43-0.4788-1.06 n D 
= 0.52, K = — = — -= 0.59 Pa s 


(600-1.703)° 5 


. , 4-6000/(1000-60) . 

v=q / A = ---r- - = \23?>m/s 

n • 0.102 2 


c) Always check the Reynolds number: 

pv d _1100 *12.3*0.102 


N =- 

iV Re 




\v 3/1 + 1" 

o— •- 

7 d An 


0.024 

Kd 
8v 


= 57 500 —> turbulent 


r 8-12.33 3-0.52 + 1 
v 0.012 ' 4-0.52 j 


0.59-0.102 

8-12.33 


= 0.024 Pas 


Use therefore a turbulent pressure loss model: 

kp = aN^-A--\pv 2 
a 2 

d) High YP is caused by slow flocculation/gelling. At long still stand (10 min or more) the gel 
strength will approach the YP. 
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Rheological models 


2.5 All models 

a) n pl = (55-40) / 511 = 0.0294 Pas = 29 cP 

x y = x - y • [x pl = 55 - 1022 • 0.0294 = 25.4 Pa 
n = 3.32 • log 55/40 = 0.459 
K = x/y = 55 / 1022 0459 = 2.28 Pas 11 

b) H & B: 1. Field approach: Take T 3 as T y , and n and K from Power law 


2. Standard: 3 points, iteration, select the 3 data points closest to actual shear 
rate, normally the three upper points 

3. Nonlinear regression: 6 data points 


c) Compare the two laminar friction expressions N? Newton ~ Bingham* an d solve for 


A Newton 1 f^eff 


d) Increased PV is an indicator of more suspended, inert fine particles (more mechanical 
friction). Increased YP is an indicator of more surface interaction between surface active 
colloidal particles. 
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Rheological models 


2.6 All models. Regression 

a) Standard 2 data points: 

Bingham: p p] = (25.4-17.2) / 511 = 0.016 Pas = 16 cP 

To = X 300 - 7300 • J-tpi = 17.2-0.016-511 = 9 Pa 

Power-law: n = log log yjy 1 = 3.32 • log 25 / 17 = 0.56 

K = x Jy* = 53 • 0.4788 / 1022 °- 56 = 0.52 Pas n 

6UU 1 oUU 

b) Regression results for Bingham (straight line) and Power law (curved line) models are presented 
in Figure 2-6: 




ou 


zU 


1U 

n 


U 

( 

i i i 

3 500 1000 1500 




Figure 2-6: Regression result as a function of# of viscometer data points. (2, 4, 6, and 8 data points). 


A summary of the results of # of data points: 


Data points 

Bingham model 

Power Law model 


PV 

YP 

R 2 

n 

K 

R 2 

2 

0.0160 

9.0 

- 

0.56 

0.52 

- 

4 

0.0192 

6.3 

0.983 

0.61 

0.39 

0.999 

6 

0.0212 

4.8 

0.977 

0.58 

0.45 

0.999 

8 

0.0232 

3.4 

0.960 

0.53 

0.59 

0.997 


The 4 upper points are the most reliable data points. The two lower ones (3 and 6 RPM) are 
influenced by the time delay inflicted by the test procedure (see Task 2.2). 
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Exercises in Drilling Fluid Engineering 


Rheological models 


cl) In the H-B model the constants are found through elimination; 


T 2 =T y + K■ Yi 
T 2 =T y +K- Y2 
T 2 =Ty+K- Ys 


25.4 = Ty + K ■ 1 022 n 
17.3= T y + K- Sll n 
13.2 — z v + K ■ 340” 


Use (1) and (2): 
to find K: 

Enter (5) in (2): 

and in (3): 

Set (6) = (7): 


r V = K- 1022 n - 25.4 = K ■ 511 n - 17.2 


K = 


25.4-17.2 
1022 n —511 n 


17.2 — T v + 


8.1 


13.2 = t v + 


y 1022 n —511 n 
8.1 


511 n 


y 1022-511 


340 n 


( —T y =) 


8,1 


t v and K are now eliminated: 


1022 n -511 n 

8.1 _ 

1022 n —511 n 


511 n — 17.2 = 


8.1 


1022-511 


340 n - 13.2 


(511 n - 340 n ) = 17.3 - 13.2 = 4.1 


(511 n - 340 n ) = (1022 n - 511 n ) 

To solve equation (10) we suggest iteration. Start with n = 0,5, and we get: 

Now try n = 0.52: 


We are moving away. Try n = 0.48: 


1.976 (22.6 - 18.4) 

= 32.0 - 22.6 


8.3 

= 9.4 

-1.1 

1.976 (25.6 - 20,7) 

= 36.7 - 2536 


9.6 

= 11,1 = 

-1.5 

3.48: 



1.976 (20.0 -16,4) 

=27.8 - 20 


7.01 

= -7.8 = 

-0.7 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


Further iteration with decreasing n to 0.35 gives a deviation cloe to 0 n = 0.35. When n is 
found, the other constants can be found easily. Now K = 3.4, resulting in Ty = - 13 Pa. 

c2) Non-linear regression. Please refer to textbook Chapter 10.6 of how to step-wise go about in 
Excel to determine the constants. 
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Rheological models 


2.7 All models. 

The rheological data are plotted in Fig. 2-7. 



We see that the Bingham model overpredicts and P-L underpredicts. This observation we test below at 
10 s -1 . 
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Exercises in Drilling Fluid Engineering 


Rheological models 


Bingham Field procedure: From the two upper viscometer data points, we obtain: 

Ppi = ^600 — ^ 3 oo = (52.2 — 35.5) ■ = 15.8 cP 

x o = 0300 - Mpi = ^ —15.8 = 17.7 lb/100ft 2 = 17.7 ■ 1.06 • 0.4788 = 9 Pa 

Bingham standard Si-procedure: 

25-17 

ju = -= 0.0157 Pas 

511 

x 0 = 25 - 0.0157 • 1022 = 8.9 Pa 

Power-law, two upper data points: 

ti — 3.32• log 25/17 = 0.556 


Herschel-Bulkley is determined through the field method: iteration and elimination, based on the three 
upper readings; 

t = r 0 + K ■ y n = 2 + 0.52 ■ y 0 - 556 


Now we test the models at y = 10 s 1 

t pl = 0.52 ■ 10°- 556 = 1.9 Pa 

t bh = 8.9 + 0.0157 ■ 10 = 9.1 Pa 

t hb = 2 + 0.52 -y 0 - 556 = 3.1 Pa 

Since the solution is 3.0 Pa the H-B is closest. The Bingham model is relatively far away, as expected, 
since theoretical model is based only on the two upper readings, which is far from the tested shear rate. 
Therefore it is, as stated in the text book, recommended to apply the data points closest to the actual 
shear rate of the model. Since these specifics are not always known ahead of time, the next best is to use 
all the data points (regression), not only the two upper ones or a few selected ones. 
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3 Drilling fluid dynamics 

3.1 Velocity profile. Continuity equation 

Define cross sectional flow area and discretize in the z-direction and then integrate velocity over the A. 


A- z-b 


dA = b-dz 


Average velocity: 


= —J az(z 0 - z)b • dz 


A 

ab 


A {) u 


_ a 


Z 0 2 1 3 

-Z-Z 

2 3 


a 

v = — 

z r\ 


f 3 


z 3 ^ 
Z 0 

3 


— CLZf\ 


1 1 


-= az a 


2 3 


v is found when the derivative of v (z)’ = 0 

max v 7 

dv _d[azz 0 -az 2 ) 


dz 


dz 


■ = az ( , - 2 az - 0 


2 az — az ( , -> max at z = — 


z o 

v = a — 

max ^ 


_ OZq 


Find v when: 


v = 8cm/s and z 0 = 4 cm 


- 2 1 
v — azg ■ - 


\ 

6 


Finally: 


8 = a ■ 16 ■ 3 
6 

8-6 12 0 , 
a = — = — = 3 m/s 

16 4 ' 

Vmax = ^~= 12 Cm/S 
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Drilling fluid dynamics 


3.2 Velocity profile. Momentum flux 


a) To find the pressure distribution: We start with Navier-Stoke in the z-direction, cylindrical 
coordinates: 


dv 1 dp u 

v — =-—+ — 

dz p dz p 


1 d(dv\ 1 d 2 v d 2 v 


r dr { dr J r 2 dO 2 dz 2 


+ g z 


Under stationary laminar flow the equation reduces to: 

0 = --^ + ^ 

p dz p 


--I 


r dr ' 

V dr)_ 


dp _ 1 d f dv 

dz ^ r dr i dr 


To replace the latter part of the equation with manageable parameters we go through 3 steps: 
Step 1: Define dv I dr froryiv^r): v = v m 


f 2 A 
1 - —y 

R 2 ; 


,v# =v r =0 


Step 2: — = v - — 

dr ^ R\ 

dv d 

Step 3: Now we derive the part r ■ — with respect to r: — r • v m 


2 r 


V max - 2 2r 

R 2 
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Drilling fluid dynamics 


The pressure distribution or pressure gradient now becomes: 


d P = .. 1 f v max; 4r 
dz r { R 2 


Total pressure loss along z: 

•4 W 


p( z ) 

\ d p=-\‘ 

p(o) o 


R 1 


dz 


^ = -\ 0 ~^ T 1 \ = 4Mv mXi z/R 2 


The shear stress is found as: r = - u • — = u v 

• • max 

or 

Maximum of x is found when r = R, i.e. at the wall as always. Navier Stoke deals with 
incompressible fluids; density does not change. 

r - 

bl) Wall shear stress: The wall shear follows from the definition of a Newtonian fluids: The x and 
r direction must be replaced by y and z respectively. 



Tw ~ T wxywall -i u ( d y+ Q*) \y=±h ~ H dy [( ( 2fi ) f 1 ^2)] \y=±h 

dv. 


T w — 




dy 


, dy 

(since v y=0 ) 



.!}ll 

' 2 y\ 

\ dz) 

' 2^\ 

K h 2 J 


-h) — 2h — 

dx 


v is reached when y = 0: 

max 1 

_ dp h 2 dp _ 2p 

Vmax ~ ^ Yx~ Vmax ' h * 


b2) Average velocity: 

It is defined as v = q/A, where q = [ vdA over the cross section. 


For our particular distribution: A = by^ dA = b dy, 



In plane Poiseulle flow between parallel plates, the average velocity is two-thirds of the 
maximum value. 
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3.3 Flow profile 

a) The equation represents the conservation of the momentum for stationary, incompressible, 
curl free, laminar pipe flow, also called the Navier Stoke equation. 

b) From momentum equation and stated condition; 

0 = -£ + «[ 1 / r -a/a> 0 ] 


From boundary condition, rearranging/integrating twice, v(r) is found. 


Alternatively, start with the shear stress vs. friction pressure for laminar pipe flow; or, one step 
later, the general shear stress equation: 


T _^_ fi H r ) 

2 dx dr 

dv{r\ = - — -r^--dr 
' ' 2 [i dx 


1 dp r 1 
2p dx 2 1 


y)=-—•—— + c , 


v(r\ = 0,C= — -R iC ^ 

' ' 4/i dx 


y(r).JL±(jp-r 

' ’ 4/i dx ' 


Dp/dx is the local pressure gradient. Since it is constant, the integrated pressure gradient is 
equal to the total pressure loss over the total length; Ap/L. 

c) Average velocity is found by taking th areal average of the integrated velocity v(r): 

- If , , 1 dp / dx R c/ 2 i \ r\ , R 2 dp / dx 

v = — \ vdA = - - - \(R 2 -r 2 )27rrdr = --- 

A J ttR 2 4p ’ 8 ju 

To compare v and v we need fist to express v ; 

■t max -t max 

v mal = v(r = 0 

By comparing the two we find an expression of axial dispersion: 

Vmax _ 1 -dp /dx-R 2 - jU _ ^ 

v 4jU-R 2 /8-dp/dx 

From this evaluation we can conclude that axial dispersion is independent of the factors shown 
above and that it is always positive; the max velocity is twice the average velocity. 
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d) The wall shear stress is found in two steps: 


Stepl: 
Step 2: 


dp/dx = 0.09 MPa/1 000 m = 90 Pa/m 
r w =S.i = °-B. 90 = 2.25 Pa 

w dx 2 2 



dp/dx 

V- 


0 . 05 2 _ 90 

8 0.0538 


0.52 m/s 


3.4a Pressure loss vs. rheology 

Average velocity is u o 


u = -j- Ju(r)dA = u 0 


A = nr 2 -» dA = 2prdr 


1 


7rR 


2 max 


1/4 


J(i-r| - 2 **-=2-^/0* 


Substitution is the solution: 




rdx = 


y = 



(l-y)R 


r/R vs. Y: 0 vs. 1, 1 vs. 0 


dr. = - dy • R 
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y m 


= |y m -(l-y)R(-dyR) = -R 2 |y m (l-y)dy 
R 2 jy m dy + R 2 Jy m+1 dy 


m+l , m+2 \ 


= r 2 ~y— + y 


m + l m + 2 


= R 


1 - 


R 


■ + 


1- 


R 


x m+2 ^ 


= R 2 


R 


m + l 


0 0 
+ • 


m + 2 


1 


m + l m + 2 
1 


' 1 1 

+ ■ 


m + l m + 2 


J J 


R 2 


R 2 


vm + 1 m + 2j (m + 2)(m + l) ^ _ _5 

4'4 


y 1/4 = 0.355 R 2 


2U„ • 0.355 R 


u = ■ 


R 


= 0.71 U m 


^=M- 


1--^ 

du du -U m V R j 


dr dr R ( r ^ /4 

1 


r —» R, t —^ go 


V Rj 

3.4b Pressure loss vs. rheology 

To find pressure loss we first need to check the Reynolds number: 


(H , = 0.016, t o = 9 Pa, n = 0.56, K = 0.52) 


a) N Re = pvd/n e ff 


v — - — -tt 2 — — 2.12 m 3 /s 

A -0.1 2 ‘ 

4 


Np„=- 


d n -v 2 n -p 


0 . 1 °- 56 - 2 . 12 1 ’ 44 -1000 


K v [(3n+l)/4n] n -8 n-1 Q 52 . p-o.se+i y^ 0f4 
V 4-0.56 J 

need to assume laminar flow 


■ = 3 500 => Turbulent, but in this exercise we 


For the Newtonian model we determine the viscosity at the estimated shear rate y = 8v/d 
8.2.12 /0.1 = 169.6 s' 1 . Here the viscosity is: 

(i = x / y = 8.6 (from Task 2.6) / 169.6 = 0.0507 Pas 
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Laminar pressure loss equations give these results: 

Newton: Ap N = 32 p L v / d 2 = 32 .0.05 07 .1000.2.12 / 0.10 = 3.4 .10 5 Pa 

Power-law: Ap pL = 4K(8v/d . (3n+l)/4n) n . 1/d . L 

= 4. 0.52 (8. 2.12 /0.1 . 2.68/2.24) 056 . 1/0.1 . 1 000 = 4.1. 10 5 Pa 

Bingham: Ap BH = 32.0.016 . 1 000.2.12 /0.1 2 + 16.1 000.9 / 3.0.1 =5.9.10 5 Pa 

We observe that Bingham “aims” high at low shear rates, and will always produce higher results 
than the Power-law model at low shear rates. 

b) Calculated shear rate for all models: 

v = 8v/d = 8 . 2.12 / 0.1 = 170 s' 1 

• Newt 

y pL = 8v/d ,(3n + l)/4n = 170 . (3 . 0.56 + 1 )/4.0.56 = 203 s - 1 

Ybh = 8v/d +T g / (3/u pl ) = 8-2.12 / 0.1 + 9 / (3-0.016) =170 + 188 = 358 s' 1 

We read from Figure 3-4 for all the models; 

x N = 8.5 Pa at 170 s' 1 : 
x pL = 10 Pa at 203 s' 1 : 
t bh = 14 Pa at 358 s' 1 : 

And the resulting pressure loss becomes: 

Ap N = 8.5 -1000 /0.1 4= 3.4 • 10 5 Pa 

Ap p L = 10 -1000/0.1 4= 4.0 ' 10 5 Pa 

Ap BH = 14 -1000/0.1 4= 5.6 • 10 5 Pa 
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Figure 3-4: Graphical presentation of viscometer data and reading of shear stress for individual shear rates. 


Model 

Y (s' 1 ) 

Read from 

Fie. 3-4 (Pa) 

Estimated from 
universal, ean (bar) 

Estimated from 

ordinarv model (bar) 

Newton 

170 

8.5 

3.4 

3.4 

Power-law 

203 

10 

4.0 

4.3 

Bingham 

358 

14 

5.2 

5.9 


One conclusion is that the two non-Newtonian ordinary and universal models give similar 
results. But the different models give different results. The best fit is as expected the P-L. 
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c) Force Balance: 

r ■ 2 7r ■ AL = Ap ■ 7r r —► r w = ^ 7? ■ ^ 

And solving Ap Newton = Ap Bingham for p eff gives the answer. 

d) A uniform pipe has an entrance length. Increasing length with increasing N Re 
Undeveloped flow at the entrance is turbulent and has higher shear until boundary layer 
expands and meets in the middle. 


3.5 Pressure loss vs. Rheology 


RPM 

f 

0 


T 

T 

rpm 

s" 1 

- 


lb/100 ft 2 

Pa 

600 

1022 

66.0 


70 

33.5 

300 

511 

47.2 


50 

23.9 

100 

170 

25.5 


27 

12.9 

6 

10 

9.4 


10 

4.8 


Mpi — 06oo — 03 oo = 66 — 47.2 — 18.8 cP — 0.019 Pas 

t 0 = 0 3OO - Hpi = 47.2 - 18.8 = 28.4 ^/t 2 = 28.4 ■ 1.06 ■ 0.4788 = 14.4 Pa 

At 33.5-23.9 00 n 

Mnz = t— = -= 0.0188 Pas 

Ay 511 

t 0 = t 300 ~ P v i ' Y — 33.5 — 0.0188 ■ 1022 = 14.3 Pa 



The Newton model can also be based on the two upper readings like the non-Newtonian. Use 
then the average value. 
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a) At this low flow rate we assume laminar flow, but still have to check; 




4 800 

n 0.109 2 ■ 1000-60 


1.425 m/s 


N re =^ = 1100 ■ 1.425 ■ 0.109/ 0.2 = 855 

V-eff 


Above we need : ii eff = \i vi + = 0.0188 + 143 °' 109 = 0.2 Pas 

r-e// r-pi ^ 6-1.425 

_ 32 \l v i-L'V 16Lr 0 _ 32 ■ 0.0188 ■ 1.425-1000 16 ■ 1000 ■ 14.3 _ n c 

bPvipe ~ ~ Slop + 3-0.109 “ /± 


b) To compare the two we need to compute Newtonian pressure loss at the correct shear rate. 


Qv 8 ■ 1.425 


YNewt d 0.109 


= 105 s' 


sv | To I 14 - 3 nr r -1 

Vdu — —— — 105 H — 356 s 

rtSM d 3 p p i 3-0.019 


Read from Figure 3-5 that x = calO and 25 Pa at the two shear rates 105 and 356 s 1 respectively. 
&Puniv,N = 4r w ■ ^ = 4 ■ 10 ~ = 3.7 ■ 10 5 Pa 
kPuniv,BH = 4- 25-^ = 9.2 -10 5 Pa 


Bingham model gives the best answers. This is natural since the mud has a typical Binghamian 
characteristic at laminar flow. We could also have compared the two by using the normal pressure loss 
models. We must then first determine the specific viscosity, just like for the universal. See also Exercise 3.4 
in this regard. 


3.6 Pressure loss. Power law 

a) Power law: t = k • y n , 

From force balance we derive: 


dv 

Y= ~Jr 


x • 2jrrL= (pi - p 2 ) tu 


-» 


A p • r 
2 L 


Using Power law model, we obtain: 


2 L 


y J^. r -n 

r 12LK J 
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Knowing that y = dv/dr we can find v(r) from: 


|- dv= l 


R r . -i- 1 

Ap 


v (r) 

v(r) = 


2LK 


r n dr 


Ap 

2LK 


n fl 


n + 1 


n +1 n +1 

R n — r n 


We need to integrate over the cross sectional area: 


q — v 7i r 


dq—v{r ) n2r dr 


-j 


= n • 2 


' Ap 1 

n n 

n +1 2n+l 

Rn _ r n 

7t • 2 • rdr = n • 2 

Ap 

n n 

2LK 

n + 1 



2LK_ 

n +1 


77+1 r) 2 

„— R n 

R n - 

2 3« + l 


• R n 


r ^ i 

i 

« 77 

i 3«+l 3/7+1 

1 p+r « . n+r 

= 7T -2 

r ^ i 

i 

” 77 

77 + 1 2/7 

_2LK_ 

77 + 1 

7 V 7 V 

2 3/7 +1 

_2LK_ 

77 + 1 

2(3n + l) 2(3n + 1) 


R 


With this expression of q we set up this relationship: 


q = n R 2 


n R 2 • v =7i 


4 p 

2LK 


• n • 


3/7 + 1 
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and solved finally for Ap: 

L 
d 


Ap = 4K 


8 v 3/2 + 1 
d 4/2 


b) Enter effective (Newtonian) viscosity for Power law fluids into the N Re expression: 


Veff,PL = ( 7 )" 1 ' K P where K p = K ■ (^)* 

pr _ p-vd _ d n ~ 1+1 p _ d n -v 2 ~ n -p _ „ 

NRe ' KpS 71 - 1 ~ K v *n-i ~ 7V Re,pipe 


c) 


Power law fluids: 


n _P _ P . r J _ J >.n . (_S_\ 2 11 - r _ 

Ka- 12”- 1 K a ( 12 »- 1 ) ^ 0 l) \^d 2 -dfj {d 2 0 -dff~ n 

(d 0 -di) n _ (do-dQ 2 ^- 1 ) 

(d 0 -d; ) 2_n (d 0 +di) 2 - n Wo +d i) 2 -" 


Newtonian fluids: 


at _ Pi v 'dh _ c i^Q (dp dj) _ Pi~ 4 • g __ ^0 _ c 

Re,Newt ^ C 2 -n-d h (d-di ) 2 p (do-dfi^d • d{) d 0 -di 


N d increases for Newtonian fluids when d. increases, while for Power law fluids it 

Re 1 

is reversed. 


The physical explanation for this behavior is the forced streamlining of the flow in narrow 
gaps. Out in the ocean, on the other hand, ocean water turns turbulent at very low shear rates. 


3.7 Pressure loss. Turbulent flow. Energy equation 

Average velocity: 


_ q 0.2 m 3 / s , . . 

v = -- = - — = 6.4 m / s 


7rR Z 


71 


(0.1 m) 


N Re = — V ~~— q qqqq^ = 128 000 -> turbulent (Kinematic viscosity; d = fi/ p) 


From Chapter 4.3 in textbook we find 8 = 0.26 mm for cast-iron pipe. Then 
s 0.26 mm 


d 200 mm 


- = 0.0013 


Enter the Moody chart (Figure 4-12 in the Textbook) on the right hand side at e/d = 0.0013 (need to 
interpolate), and move to the left to intersect with N Re = 128 000. Read: f = 0.0225. The head loss becomes: 





= (0.0225) 


500 m 
0.2 m 


(6.4 m / s) 2 
2(9.81 m/s 2 ) 


= 117 m 
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3.8 Pressure loss vs. flow rate 

From exercise 2-7 we copy these data and present the flow curve: 
n = 0.56, K = 0.52 Pa s 056 



Shear stress (1/s) 


Figure 3-8: Flow curve. 


The pressure loss equations: 


Ap lam = 4 K (— ■ —)* - ■ AL = 4 ■ 0.52 (— ■ —)°' 56 ± ■ 1000 = 267 553 

rLam \d 4 K ) d VO.l 2.24/ 0.1 


,0.56 


AT 1 

^ = aW-»— --pv 2 =1.619-10 6 W E 

d h 2 


-0.214 2 


a = (log n + 3.93) / 50 = (log 0.56 + 3.93)/ 50 = 0.0736, b = 0.286 
b = (1.75 + log n ) 17 = (1.75 + log 0.56)/ 7 = 0.286 


V = 


g_ 

A 


q 

7T /4 • 0.T 


= 127.4 q 


M Re 


p-vd _ d n -v 2 n -p 
R-eff Kp-8 n 1 


0 


0 . 1°- 56 ■ v 1 - 44 -1100 

.56 • ° ^ 6 1 8 _ 0.44 

V 4 0.56 J 


1 206 ■ V 1 - 44 
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Select three data points in the laminar and three in higher flow/turbulent regime and plot the result 
(Figure 3-8): 
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Figure 3-8: Laminar and turbulent pressure loss. 

3.9 Pressure loss. Field data 

Check first the solution for 5 spm. Then repeat for the remaining speeds. Here a check is made only 
for the 5 spm case. Here we think the flow is laminar. Pressure reading at 5 spm =13 bars. To check 
theoretical Ap at 5 spm we go through these steps: 


Step 1: Rheology: n = 0.50, K= 1.63 Pas' n 


Step 2: Check flow regime: N Re 


pvd 

Veff 


Peff,l 


/ 8v 

3 n + l\ n 

Kd 

/8 ■ 0.15 

3 ■ 0.5 + l^ 0,5 1.63 

■0.105 

Vd" 

4n ) 

8v 

V 0.105 

CD 

LO 

o 

0.15 


= 0.54 Pas 


Peff.Ann 


/ 12v 2 n + l\ n Kd 

/12-0.02 2 ■ 0.5 + 1\ 

05 1.63-0.188 

\d 0 — d { 3n / 12i; 

l 0.188 3-0.5 J 

1 12-0.02 


= 1.47 Pas 


v = - A , d = 4.127 -0.0254 = 0.105 m 
A 

q = 15.29 ■ 1(T 3 ■ 5 = 0.0013 m 2 /s 
A 0 p = n / 4 ■ 0.105 2 = 0.0087 m/s 
A Ann = (12.40 2 - 5 2 ) ■ 0.0254 2 = 0.065 m 2 
v d = q /. = 0.0013/0.0087 = 0.15 m/s 

F / /i op 

v Ann = q / A = 0.0013/0.065 = 0.02 m/s 

' 71 nn 

do - di = (12.4 - 5 ) 0.0254 = 0.188 m 


p ■ d ■ v 

NRe,DP = 

Peff 


1120-0.15-0.105 
054 


= 32 (since it is so low it is assumed laminar also in the annulus ) 
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Step 3: Find total pressure drop in flow system: Ap = Ap dp + Ap ann 

L 4500 

A Vd P = 32vn eff — = 32 ■ 0.15 ■ 0.54 ■ - iq52 = 10.6 ■ 10 5 Pa 
L 4500 

A Vann = 48 Vfi eff — = 32 ■ 0.02 ■ 0.147 ■ - = 1.2 ■ 10 5 Pa 

A Ptot = A Pdp + A Vann = (10.6 + 1.2)10 5 = 11.8 ■ 10 5 Pa 

The task was to evaluate and compare the estimated pressure with the recorded pump pressure. As 
exemplified at 5 spm flow, the theory compared to reality, 11.8 vs. 13, differed by -9%. One possible 
explanation for this difference is that the theoretical Power law underestimates the real shear rate. Other 
error sources exist; selecting the correct rheology model; fluid is very Binghamian. When does flow turn 
from laminar to turbulence? In real wells it turns turbulent earlier than theoretically predicted, due to 
roughness and uneven flow path. Turbulent pressure loss is highly empirical. 

3.10 Pressure loss. Effects of rotation 

During rotation of the pipe, external energy is added to the fluid (from the top drive). If the pipe is 
exposed to strict, controlled rotation the mud will be stirred (like in a viscometer). The mud is non- 
Newtonian and shear-thinning (n < 1.0) during laminar flow. For laminar flow at SPM = 5, the friction 
loss will be reduced from 13 at pure flow to 8 bars at added rotation. 

To check this fact we will try to estimate the reduction in effective viscosity. Assume the fluid behaves 
according to the Power-Law. 

Step 1: Find v rot = peripherical velocity of pipe. One rotation corresponds to the length 2jtr. The outer 
wall velocity is zero. Divide by 2 to obtain the average of the two extreme velocities: 

/ min \ 100 1 

v = [2nr- RPM ■ 1 ——s + 0 /2 = 2n ■ 5/2 ■ 0.0254 ■-—■- = 0.33 m/s 
\ 60 / 60 2 

Step 2: Translate the rotational velocity variation to shear rate: 

12v 2n + l ( 12-0.33 2 ■ 0.5 +1\ 

y =-= [-1 = 30 S -1 

Y d 0 — d t 3n V13.4- 5) -0.0254 3-0.5 / 

Step 3: Calculate the shear rate caused by pure axial flow at 5 SPM: v ann at 5 spm was found to be 0.02 
m/s in the solution of Task 3.9 above. 

12v 2n + 1 _ / 12-0.02 2 ■ 0.5 + 1\ _ 

7 ~ d 0 ' 3n “ V13.4- 5) -0.0254 ’ 3-0.5 J ~ L5 s 
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Step 4: Effective viscosity: In Figure 3-10 the shear rate has been increased from 1.5 to additionally 30 
s 1 (the two were assumed to be additive). We compare the effect of rotation in terms of changes 
in the Newtonian viscosity: Use the flow curve in Figure 3-10. Here we are left to assume the 
shear stress at these low shear rate values. From the rheological data in Task 2-2 we select the 
smallest recordable one, 3 Pa. 

u =t/v = 1/1.5 = 0.67 Pas 

k 5 spm * 

li tf . tinn = x / v = 8 / (1.5 + 30) = 0.25 Pas 

r 5 spm+rotation at 100 rpm * v J 

We see that the effective viscosity has reduced to less than half of its original at that flow rate. 
Note also that the shear stress seen from the pump (SPP) is not increased, it is “for free”, only 
the effective viscosity is decreased, as shown in Figure 3-10. 
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Shear rate (s _1 ) 

Figure 3-10: Effect of rotations on the rheology results in new ECD for instance at the low shear stress of 10 Pa. 


If the pipe is exposed to vigorous rotation and possibly vibrations, which often is the case in a real well, 
then additional turbulence will introduce higher pressure loss, not a lower one. However, at the highest 
pump speed (100 spm) flow is already turbulent and the complex motion of the drill string will only 
negligibly influence the pressure. 

3.11 Pressure loss. Bit nozzle. OFU 

The exercise information is summarized below: 



a) In order to estimate the annular pressure drop, the mud velocity, the effective viscosity and the 
Reynolds number are needed, valid for the BHA-annulus, while drilling at a depth of 2500 m: 

Gal 0.1311ft 3 / gal 

■ - ™ 2 ,4 „/, 


V=«: 

A 


min 60 s / min 


= 0.40856-^ = 
dl 


-d 


in • 


1 


\2 


12 in! ft 


2.448 d\ 2.448 (l2.25 2 - 6.25 2 ) 
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First check the flow regime before estimating Ap: 


^Re = 


__ lSlpv{d 2 -d^) _ 757• 12.9• 2.94(12.25-6.25) _ 


= 945 => not turbulent 






5 r-d b 5-15 (12.25 —6.25) 

ll„ = u p , +—-—- = 35 + - -- - = 119 cP 

eff PL v 2.94 — 


Ap Mpr} 


A l dc \tm(d 2 -d,f 200 (d 2 -di) 


— -= 0.0154 psi/ft => 0.0154- =5.05 psi 

0.3048 


Eqv. p in 2500 MD: 

Ap, 


Pek V =P + 


sm — = \2.9 + 5,05 =13.15 PPG 


0.052 -L 


0.052- 


2500 


0.3048 

b) Bit pressure drop in OFU and in Si-units: 

8.311-10 5 -12.9-800' 


APbit- 


0.95 • 


5---I — 


V 


4 132 


2\ 


1 346 psi 


A p bit = 1.11 -1550 


800 • 


3.78533 
60 1000 


n ( 14 


4 l 32 


• 0.0254 


= 9.318 MPa = 1369 psi 


Conclusion: Choice of units does not make any difference for this formula. 


c) Total pressure loss in the system: 

= 17.87 ft Is 

r 2 _ J I 


DP: 

= _ 800 

V_ 2.448-4.276 : 

DC: 

800 

V _ 2.448-3.1 2 ’ 


: 34.0 ftp 

Turbulent or laminar in largest area? 

Ar _ 928/OvJ _ 928-12.9-17.87-4.276 

■™ Re — — TT? 

Peff 35 + 

Pressure losses, turbulent: 


= 26136 ^ turbulent by far 


0.75 .. 0.25-1.75 


DP: A Plurh / Al = 


P P P i 


12.9 0 ' 75 -17.87 1 ° -85 


1.75 oc0.25 


1800<i 125 1800- 4.276 125 

Total pressure loss in the circulating system (not all details shown): 


■ = 0.29 psi / ft 


A P tot = A PdP + A Pdc + A Pbi, + A P annl + A Pan„2 

Ap tot = 2 320 + 300 + 1 346 + 5 + 100 = 4 071 psi 
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3.12 Swab pressure. Clinging factor 

a) Make a detailed drawing. The clinging fluid is defined by R in Figure 3-12, causing a certain 
volume to accompany the string upwards. 

b) Identify an expression of the velocity profile with the correct boundary condition; 

v ( r Uo = ° 

V ( r Ucling = 0 

y ( r Ui = V p 

From the velocity profile we can obtain the cling volume pr. unit time: 

qcllne = j v(r)2irrdr 

Rcling 

c) Mass balance: q = v . • A + q ,. . An equivalent fluid mass has to move down 

' T-up pipe pipe T-ciing ^ 

d) Cling factor = q .. / q, 

7 o J-cling T-down 
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Figure 3-12: Flow streaming down the annulus when pulling a closed-end drill string. 
The frictional pressure loss causes a pressure deficit below the string's bottom end. 


3.13 Swab pressure model 

The V dm will be added to the fluid volume displaced by the drill string. 


^DS ^DS ‘ V DS’ ^displaced ^DS^flow 

v = V _ V _ V 

total pump displaced cling 

Step 1. Understand the physics and forces in a Drawing 
Step 2. Make a control box and entering exiting forces. 

Ap-Jtr 2 = x-2jtrAL 
Ap/AL -r = 2 t 

Step 3. Differentiate: Ap/AL • dr = 2 dx 

Step 4. Integrate over the control volume (along r) 
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4 Hydraulic program 

4.1 Mud pump issues 

a) Characteristic diagrams of mud pumps are shown in Figure 4-1: 




Figure 4-1 : Maximum pump pressure for each liner (left) and maximum effect for the same (right) 

Power or effect = E = p • q, is expressed in Horse Power (HP) or in watt (W). As shown in 
Figure 4-1 (to the right) the maximum effect is constant for all liners. 

Volumetric efficiency = 0.96 for new pumps. When it has reduced to 0.93 the pump has to be 
upgraded (replace valves and / or piston / liner). 

b) Parallel —> higher q, used for upper wellbore sections. Series —> higher p (only for very special 
application, not done for mud pumps). 

c) Centrifugal pumps can deliver high flow rate but limited pressure (max 3 bars) since there is 
an open slot between statot and rotor. 

d) Knocking is a result of high acceleration. At the highest acceleration the water is not able to 
fully follow the high accelleration, and the following takes place: -> vacuum effect -> vapor 
pressure is surpassed and water starts boiling (vapor bubbles are formed on the steel surface) 
-> implosion of vapor bubbles during retardation of piston movement -> repeated shock vaves 
weakenes the steel. 
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e) Find first A p d = A p p - A p bit 

Bit pressure loss can be determined since nozzle sizes are known: 
p bit =\.\\^pv 2= 1-11 1 ~' 1300 ■ 100 2 = 72.15 ■ 10 5 Pa 

V C \! Anozzle 


Anozzies= (Tt/4)d c 2 = tt/ 4 ■ 0.02 06 2 = 3.3 3 3 3 • 10" 4 m 2 

d e = Jd i 2 + J 2 2 +J 3 2 =73 -d x =73—0.0254 = 0.0206 m 

= 2 500 /(I 000 • 60) = 0.04167 m 3 /s = q/A = 0.04167 / 3.33-10 4 = 125.2 m/s 
q 2 = 670 /(I 000 • 60) = 0.011 m 3 /s v 2 = q 2 /A = 0.017 / 3.33-10 4 = 33.5 m/s 

Resulting in these values: 


Parameter 

q 

Pp 

A p bi , 

AP d 

Unit 

Ipm 

bar 

bar 

bar 

Data 

2500 

230 

72 

157.9 


670 

33 

8.1 

24.9 


K x and m are determined from: Ap d = KjD q m . Selecting the operator In or log does not make 
any diffidence here. 

m = In (Apd 1 /Apd 2 ) / In (q/q 2 ) = In (157.9 / 24.9) / In (2 500 / 670) = 1.6 
K x = Ap d /q m • D = 157.9 • 10 5 / (0.04167 1 - 69 - 2 500) = 2 • 10 6 

4.2 Nozzle selection. Section wise 

a) Existing parasitic pressure loss and its parameters (assume 3 identical nozzles): 

d e = Vn ■ d,. = 73 g ■ 0.0254 = 0.0206 m 

= * =_= 100.2 mis 

A 1000-60- —0.0206 2 0 ' 000333 

4 

Vnozzle,2 87.7 771/S 

Vmti ~ 1-11 ■ 1/2 • 1400 ■ 100 2 = 78 ■ 10 5 Pa / 59.8 ■ 10 5 Pa (represents pump rate 1/rate 2) 
Apioss = 186 - 121.8 = 108.2 bar 
Apioss = 120 - 59.8 = 60.2 
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A P,oss =K l D-q‘ 


m _ logCAP^/AP,^) _ lofe< 108.2/60.2) . lM 

lo &qjq{) log( 2500/1750) 


1 2500-0.0333 164 

b) From existing hydraulic program parameters we obtain: 

%t = [ 2 Ppumpm / ( m+2 ) K i D ] 1/m = [2-270-10 5 ( 1.64+2)-1.16-10 6 -3 000)] m - 64 = 0.0043 °- 61 = 
0.0358 m 3 /s 

Check the flow rate at the boundary conditions (in Si-units to be able to compare): 

Qmax.p = 2 700 Ipm = 2 700/(1 000 • 60) = 0.045 m 3 /s 
q r = 1 700/(1 000 ■ 60) = 0.028 m 3 /s 
= 2 600 / 60 000 0.0433 m 3 / 5 

Since q is in between q and q , 0.0358 is the optimal solution at this depth. 

T-opt 7L r 7L ma x,ann * r 


Maastricht University 
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c) The optimal bit pressure loss is given by this equation: 

Ap bit = Pmax.p ■ m/(m + 2) = 0.458 ■ 270 ■ 10 5 ■ 1.64(3.64) = 121.6 ■ 10 5 Pa 


To find the optimal nozzle, solve pressure drop equation vs. optimal bit pressure drop 


l.ll--p(ff) = 121.6 -10 5 Pa 

2 V A ^e) 


di = 


1.11- 0.5 p-q 2 


(?) 


121.6 - 10 5 


j4 4 1 . 11 - 0 . 5-1 400 - 4 2 - 0 . 0358 2 

dl, = -----= 0.019 m 

e \ n 2 ■ 121.6 -io 5 


Nozzle diameter is not optimal at 3000 meters depth but close to. It should have been: 


, d e 0.019 

d S ,=-r i = 


j 0.019 32 „ „ . . 32 . „ „ ... . N 

dnFjj — -?=-' _ — 0.011 ■ _ — 13.9 (32nds in) 


V3 V3 U ° FU V3 0.0254 


0.0254 


To stay below maximum pump pressure we select: 


14/32", 14/32", 14/32" 

4.3 Liner selection. Section wise 

a) To optimize q opt you need criteria; maximize ROP. First find a suitable expression of ROP: 

ROP=A-(q/d nozzle r 8 


Further we need to know the boundaries: 

Pp = Apio S s + ^Pbit 

Pp = K 1 Dq m + K m ■ (£) 2 -» 

p p = K 1 Dq m + K bit -(jV) I q 2 

Pp ■ q 2 = KiDq m+2 + K m ■ 

± = {Pt ,2_ kiDqm »y^i_ 

a e K bit 

Entering the ezpression of q/d into the ROP expression, differentiate and fid maximum when 
dROP/dq = 0, resulting in an optimal expression 

2p p q-K l D(m + 2)q m+i =0 
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b) In order to cover the complete range of possible flow rates, we see from the data and the pump 
characteristics that the 6%" liner has a q max of 0.0422. Only this liner is large enough to supply 
the q max of 0.04 m 3 /s if necessary. Intuitively we would select the smallest liner possible, the 
5 V 2 ", since it has the highest available pressure. This is what we also recommend. Alternatively 
we could check what optimal flow rate each liner produced at this depth: 


tfopt5.5 - 


AP'l 

m 

2-383 -10 5 

K, ■ D(m + 2) 


1.6-10 6 -2500 (1.7 + 2)_ 


Qopt6 


Qopt6.5 


optl.5 


2 (p p l 

I 

m 

1 

2-322-10 5 "p 

K r D(m + 2) 


1.6 -10 6 - 2500 (1.7 + 2) 

2 (Pp) m 

1 

m 

2-274-10 5 1 T 

K r D(m + 2) 


1.6-10 6 -2500(1.7 + 2) 

/ 2 (pX 

1 

m 

2-220-10 5 l 1 

K r D(m + 2) 


1.6-10 6 -2500(1.7 + 2) 


= 0.0454m 3 Is 

0.0405 m 3 /s 

= 0.0369m 3 Is 

= 0.0324m 3 Is 


We observe that the optimal flow rate is outside the boundaries for 5.5, 5.75 and 6. We could 
claim that the 6.25" liner is the smallest liner with its optimm within the rang and therefore 
select that one. 


An expression of optimal bit pressure is found from the same optimal solution as above. 
Substitute p d = K x D q m : 


2 P p ~(m + 2)p d =0 
2p p -(m + 2\p p -p bit )=0 

v in 

Optimum bit pressure drop is thus: p bit = -— p p 

m + 2 

From pump data of 6" liner we read: p = 322 • 10 5 Pa 

r r r pump, max 

Pbu,opt = F ■ 322 ■ 105pa = 150 • 105 pa 


c) The depth of altering liner is found when q opt ^ of 6" reaches q max of liner 5% " which is 0.0307 m 3 /s 


2 _jpi L 

K r D(m + 2) 


0.0307 


K r D(m + 2) = 


2 (^L 

(0.0307) m 


2 [Pp\ m = (2 * 322) / [1.6 * 10 6 * (3.7 * (0.0307) L7 ] = 3 841 m 
K r (m + 2) (0.0307)" 


91 


Download free eBooks at bookboon.com 

































gence cdg . © Photononstop 


Exercises in Drilling Fluid Engineering 


Hydraulic program 


4.4 Hydraulic program. Section wise 

a) By including all the variables in the ROP equation we might succeed. 


ROP = A (q/d e ) a s 

Pp = A Ploss + A Pbit 


p p = K^Dq'" +1.11- 0.5p 
P p = K x Dq m + k bit 


\ 2 

J__ 

AO 


1 


<? 

kPj 


{n / 4) 


= 1.11.0.5. p / (jt/4) 2 


j- = {p P -K Dc i m V-~jjj 

ft 


Substitute into the ROP equation: ROP = A 


-\a % 

( Pp -k 1 D q '"J 5 .-L =(-Lf 'A-(p p —k 1 Dq m J S 


From this equation we see that for increasing depth the ROP-function will decrease. 




> Apply now 


REDEFINE YOUR FUTURE 

AXA GLOBAL GRADUATE 
PROGRAM 2015 


redefining /standards 


92 

Download free eBooks at bookboon.com 
















Exercises in Drilling Fluid Engineering 


Hydraulic program 


b) Only the case of vertical wells is presented here. The horizontal case is solved identically, only 
the flow rate has higher importance / higher upper boundaries. Within the narrow flow rate 
interval for vertical wells, the 6 l A" liner covers the maximum flow rate. First check what is the 
optimum flow rate for each of the liners, from 6 l A" and downwards at 1 950 mMD. Since the 
boundary conditions do not change with depth, intuitively we would suggest that the smallest 
liner possible should be selected, which we do; the 5V£" liner. Just for the interest we also check 
the optimal flow rate at the two depths: 


Qopt, 1950,5.5" 
Qopt, 4 000,5.5" 


= r 2-App 1 


1/m 


l(m+2)K 1 Dl 

■ 383■10 5 


-k 


2•383■10 5 


= [- L 

L(i.65 +; 


2)2 ■ 10 6 -4 000 


(1.65 + 2)2 • 10 6 ■ 1950 
-.1/1.65 

J = 0.0273 m 3 /s 


,1/1.65 

J = 0.0418 m 3 /s 


We observe that when drilling at 4 000 m the hydraulic program is inside the theoretical limits, 
but not while drilling at 1 950 m. Our suggestion is therefore to use the q max , 0.028, until a 
depth at which q opt is reached: 


D = 


2 ■ A Pp 

(m+ 2 ) K x q m 


_2.383 ■ 1 0 5 _ 

(1.65 + 2) 2 ■ 10 6 - 0.028 1 - 65 


= 3 829 m 


Then use q opt till TD. If change of liner is allowed with increasing depth, then a better solution 
is to select the liner corresponding to optimum flow rate. 


4.5 Optimal parameters with BHHP. OFU. Section wise 

Optimum parameters for this criterion can be derived as follows: 

ROP = C ■ BHHP = ^ ■ A p bit ■q = c 1 (p p - A p loss ) ■ q 
ROP = c x (p p ■ q — K-tDq™ ■ q) 


Only the derivation of the core produces a sensible answer. 

= p p -(m + l)K 1 Dq m = 0 

Ap p , op t = (m + l)K 1 Dq m => q opt = [ ^.^J ^ 

We observe also that optimal parasitic and bit pressure drop can be derived from the equation above: 
A Ploss ~ Pp ^Pi bit,opt ~ m+1 Pp 
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To estimate q opt , K X D and m must be determined from: A p loss = K x Dq m 


Input 

Output 

q (GPM) 

P p (psi) 

A Pbi,(P si ) 

Ap d (psi) 

500 

3 000 

2 097.0 

903.0 

250 

800 

524.3 

275.7 


_ In (hpiossi/hpiossz) _ ln(903/275.7) _ 

\n(qi/q 2 ) ln(500/250) 

KiD = Apioss/q m = = 0.02168 Q psi/GPM m ) 


Now q opt can be found (note that the units are correct): 


_ [ &Pp.max ] 1/m _ \ 

Lcm+l^oJ Li 


3620 


11/1.712 


(1.71 + 1) 0.0217J 


= 627 GPM 


Including the volumetric efficiency factor, e vol , the practical q opt is: 

q opt = 627 • 0.9 = 564 GPM 


Now the bit nozzles can be estimated from Ap bitopt = Ap bit : 


m+l ' Pp> max ~ 5 ' P ' Q /Ql ' ^nozzle 

_ j8,311-10- 5 -pql pt _ |8.311-10- 5 -564 2 

A nozzle,copt ~ ^ c 2 . Apwt “ 0 .95 2 -2284 

Pbitovt = — ■ 3620 = ■ 3620 = 2284 psi 

ran,opt m+1 1.71+1 p 


= 0.33 in 2 


Finally the area can be expressed as: 


A = 



dl 


Tl-4-32 2 

ttV3 


d e = J-^j \ ^ ~ x (^2nd in) 

Assume x = 15.8 ^ 16/32", 16/32", 16/32" 
Assume x = 16.8 -> 16/32", 17/32", 17/32" 
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4.6 Liner selection. Complete well 

a) Optimal flow rate in pump range II is derived initially like in Exercise 4.4. 


ROP = K • 


r \ 

£ 

V dnozzle J 


~{p P q 2 -K x D- q ^f 


K, 


In this equation we substitute qp = E 

ROP„ = ^- t (Pmax ■ q 2 ~ Kt D ■ q m+2 p/ 4 = K'(E V ■ q - K l Dq m+2 f*U 


8ROP n _ K ] (E p - (m + 2) K x Dq m+> ) < 


dq 


, A 

( ) 4 


= 0 


Q opt a 


K r D(m + 2) 


i 

m +1 


b) Range II stops towards lower q, when q = q 5 5 „ max , while range I stops towards higher q when 

flopt,I fl5,5”,max. 


Qopt, 


(^£>(771 + 2)) 
1/m 


l/(m+l) 


Du = 1852 m 


1/1.5 


/ 2pp / 2-383-10 5 ono 

Qcmt i = 7-f — = - — = 0.028 =3 Dj = 2123 m 

V(m+2)/C 1 D/ \3.5-2.20T0 6 -D/ 1 
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c) The graphical solution speaks for itself: 



Figure 4-6: Graphical solution. 


d) From Table 4-6 in the Task section we can see that between 1 500 and 3 000 m there are two 
potential options. However, from Figure 4-6 above we already know that at 1 852 m we are 
inside pump range I, and applying the SVi" liner. We also see that the maximum flow rate of 
the 5 Vi" liner is the optimum one for most of its operational time. This has already been shown 
in task b) above. 


4.7 Liner selection. Complete well 

a) Find first what working area prevails in the boundary zone between pump are II and I, which 
are valid as indicated in Figure 4-7.1. Area II ends with the 5.75" liner when q max 55 is reached: 


Ropt.n ~ 0.028 — [- 


-il/(m+l) 


r 1.076 -10 6 I 1 / 2 - 52 

= --- = Djtt = 1457 m 

Ll.72 -10 6 - D • 3.52J 11 1 


X!D(m+2)J 


At 1 457 m depth the pump area I is entered, coming from II. From the pump chart we note 
that 5 Vi" liner has p max = 383 bar. Solving q optl eqn. for D we find where q optI starts. 


Roptj 


0.028 


paH 1/m = [ p" = = 2 900 m 

L/< r 1 Z)(7n+2)J Ll.72 -10 6 ■ D ■ 3.52J 1 ° vz 


1/1.52 
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b) In Figure 4-7.1 we see the division between the two pump areas. A preliminary solution of 
how to distinguish between them is given in Figure 4-7.2, its upper part. 



Figure 4-7.1 : Boundary area between pump area II and I 



Figure 4-7.2: Flow chart of estimating hydraulic parameters while in the boundary zone between the pump area ranges. 
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5 Wellbore challenges 

5.1 Filtration control 

a) Fluid loss can be optimized in three ways: 

1. Ensure a good filter in the wellbore wall by including a wide range of particles in the 
mud (Bentonite, Barite and cuttings are already there. It may be enough) 

2. Ensure a good filter cake by making sure the mud is dispersed 

3. Losses through the filter cake in WBM are treated by increasing the viscosity of the 
fluid phase (starch, CMC and sodium polyacrylate). OBM rarely requires such additives 

b) Increased fluid loss is detected through the filter press test. Depending on the fluid loss, its 
rheology and other test results you normally find out which of the 3 qualities in question a) 
above are the problem. See more details in textbook Chapter 2. 

c) Filtrate will lead to approximately the same invasion depth (typically 3 feet). Particles will invade 
both sands typically 1 ft, assuming the particle size distribution of the added or the existing 
particles in the mud covers at least the range of l/7 th to l/3 rd of the pore size distribution. 
That assumption is normally fulfilled through the presence of Bentonite / polymers + Barite 
+ cuttings. 
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d) Water-flow into the shale can be controlled by 

1. Increasing the viscosity of thefluid face 

2. Establishing a filter; same principles as in c) above 

5.2 Filtration control 

a) Standard fluid loss = filtrate volume in ml after 30 min 

V f =AC-4t 

API filter press A x = 45 cm 2 = 45 • 10' 4 m 2 
Flow into porous formation: 

Effective formation height (10%) h = 4000 ft • 0.3048 • 0.10 = 122 m 
Wellbore surface area A 2 = jt-d*L=15- 0.0254 • n • 122 = 149 m 2 

Assuming that an identical filter is built up on the sandstone surface, the loss during the first 
30 minutes can be estimated: 




The filter will have the same progress and same constants. 


b) Intrusion after 24 h with 100% displacement of 5.74 m 3 . Since only the pores are invaded, the 
affected volume to be invaded will be accordingly larger. 



eJJ 0.15 

Hole radius: r hole = ^ • 0.0254 = 0.19 m 
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Effected or invaded volume can be expressed by means of the radius, r. , , 

J invaded 

V eff /0 = (x-r invaded -x-r h o le 2 )-h 

invaded 2 ^ = V ff /0 + Jt • h 

R i„v=^ # /0 + ^- r,jyil{nh) 

R inv= 75.74/0.15 4- 7 1 ■ 0.59 2 - 122)/(7t ■ 122) = 0.369 m (which is round 1 ft) 

c) Obtain knowledge of the sandstones pore size distribution. Apply 

• Particles in the mud with a mean size = V 3 of pore size. 

• Colloidal particles must be dispersed 

• Base fluid viscosity can be increased. 

For OBM the particle distribution is manly taken care of by emulsified water droplets, with a 
wide droplet distribution. 

Follow up by making a filter press test every 30 min. In case filter loss increases, every bullet point 
above must be evaluated through more tests, viscosity especially and specific countermeasures 
are set in motion. 

Filter test is important mainly to obtain a thin and tight filter cake to avoid stuck. Secondly, 
while drilling in the reservoir, to minimize fm-damage. 

Filter loss is increasing because we are drilling trough contaminants like salt, chalk, etc., which are 
causing flocculation; resulting in thick filter cake and perhaps in stuck pipe. 

Solution/explanation: Add thinners/-dispersants/deflocculators. 

If massive attack; dump the returning mud. 

5.3 Cuttings concentration 

a) At the end of the horizontal section the concentration would also be close to 0.02. A bed has 
been established which is in equilibrium between erosional and depositional forces. The bed 
height is a result of the steady state process. Settling and lifting of cuttings is a function of 
flow rate, ROP, etc. The bed height can be estimated theoretically. We select not to do that, but 
rather assume that a certain % of the generated cuttings are left behind in the wellbore. E.g. 
10% of the wellbore cross sectional area is in average filled with cuttings. This means 10% of 
what is drilled out is left in the well. The initial concentration is accordingly reduced by 10%; 

c 2 = c 1 . (1-0.1) = 0.2.0.9 = 0.018 
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b) In the vertical section the concentration will increase due to the slipping of cuttings, expressed 
through the transport ratio R t . This could have been estimated through Stokes Law. We select 
here just to assume a typical value of R t , which is 0.75. We then have: 

c 3 = c 2 / R t = 0.018 / 0.75 = 0.024 

c) The cuttings removal forces are opposed by gravity and cohesive forces. Gravitational forces 
are given by Stokes law. Removal forces are given by drag and lift forces. The drag force is a 
function of particle Reynolds number and spherisity. The drag and lift forces can be developed 
info the critical lift velocity or rolling velocity. String RPM will induce a drag force onto the 
cuttings, referred to as the viscous coupling. At a given cuttings feed rate, a stationary bed 
height will form as a function of the equilibrium state between all involved input variables, 
some eroding, some building. 

d) Here the bed height will be higher due to higher cross sectional area and thus lower fluid 
velocity, assuming the entrance effects do not ruin this assumption. Settling rate will increase, 
but first after the entrance effects are passed. When BHA during tripping is shoveling much 
cuttings into a narrower wellbore it is easy to imagine that the BHA may become jammed. 



Brain power 
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e) 


At least five issues are related to poor solids control: 

• Building of solids bed —> Mechanical stuck pipe 

• Building of viscosity / high gel —> high pressure losses -> high ECD —> fracture 

• High cuttings concentration in annulus —> higher mud density in annulus —> high 
ECD —> lower ROP 

• Higher mechanical friction —> loss of WOB —> reduced ROP 

• Loss of filter control —> thicker mud cake —> differential sticking 


5.4 Cuttings 

a) Friction, drag, lifts cohesive, gravity Make drawing. Rolling when > 0, lifting when 

b) Stokes law is valid only for infinite dilution (one spheres), for small spheres (r < 0.1mm) and 
for laminar flow around the particle. 

c) High ROP, high inclination, highly viscous mud, rolling, lifting etc. 

d) Two drawings are needed: One for Drag, Lift, Cohesion, Gravity, and one for a rotating string 
with viscous coupling. 

The viscous coupling radius relative to R, increase with RPM and climb the TJ at around 120. 


5.5 Density control 

a) When weight material has to be added to the mud it is called weighted mud. To clean it, 
normally a 60 + 250 Mesh shakers + a centrifuge would be perfect. Barite is preliminary taken 
out by the centrifuge and later re-injected. 

b) By combining this information: 

'"-"■r, 


Yj m and m, = V/p. we obtain 

Pl 

Using it for the specified task: 


• AV„ = — 


Pb Pi 


AV _ 60000(1.4-1.3) _ 
B 4.3-1.4 


= 2 092 / 


However, barite is controlled by its mass, not volume: 

A m B = AV B -p B = 2.092 m 3 • 4.3 ton /m 3 =9.1 ton 


c) In the next two tasks we simply take the weighted average of all densities: 

- PZi + P£l + PjLl + A,' K 1500-1 + 1800-0.1 + 40 +m B -(m B /4300) 

P ~Y/~ V i+ V 2+ V 3+K ~ 1 + 0.1 + 40/2300 + m B /4300 

m B = 143.2 kg 
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d) = 1.55 = 


P\ V \ + Pl V 2 + PtVs + Pw' V w 


v l + v 2 + v 3 + v w 

1500-10 + 1600-20 + 1900-3 + r w -1000 
10 + 20 + 3 + F, 


= 1550 -> 


V = 2.8 m 3 


5.6 Density control 


a) Remaining volume v r = ^ V = 60 m 3 


Remove 40 m 3 and then add = v w + v b = 40 m 3 


Density of added volume must also be 1.8 kg/1 

= m 1 + w w + m, 

Padd K + V w+ V b 


m x = /?,- Fj= 1800 • 60 = 108 000 kg 

V w +V b =40 m 3 


mo = 108000 + ™^ , 
100 


m + m = 72 000 kg 

w b ® 

m = V • 1 000 

w w 

m b = V b -4 300 = (40 - VJ 4 300 
V w -1 000 + (40 - V w ) 4300 = 72 000 kg 
V = 30.3 m 3 

w 

V b = 40 - 30.3 = 9.7 m 3 

m b = 72 000 - m w = 72 000 - 30.3 • 1 000 = 41 700 kg 


Conclusion: Remove 40 m 3 mud, add 30.3 m 3 of water and 41 700 kg Barite to obtain 100 m 3 
mud of p = 1800 kg/m 3 and LSGS of 3%. 

b) Dump first V 2 = 90-(0.055-0.035)/0.055 = 40.9 m 3 

Remaining volume Vj = 49.1, Pi = 1 600, ny = 78 560 kg 

p w = 1 000, p b = 4 300 kg 

Step 1: Add 40 900 kg of water and find new density: 

p new = (78 560 + 40 900) / 90 = 1 327 kg/m 3 = p old in step 2. 


Step 2: Find out how much Barite must be added to reach the desired density of 1.7 kg/1 


v - v 

V add v 1 


Pnew Pold 

Padd Pnew 


= 90-(1.6-1.327)/(4.3-1.6) = 9.1 m 3 
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The new volume fraction of Low gravity solids, becomes: 


LSGC 2 = 0.03 * (90 +9.1)790 = 0.033 3.3% 


c) Necessary r2: r 2 gh = 410 • 10 5 —» 


410 10 5 
9.81-3000 


= 1393 kg/m 3 


Fill a tank of 10 m 3 with new mud of p 2 : 


From data: 


V 1 = volume old mud V 2 = 10 m 3 
p l = 1 200 kg/m 3 V b = volume barite 


p b = 4 300 kg/m 3 


Increase from pl to p2: 


Pi 


v b = h ? Pi = 


V, 


F 1 + V b - P b 
V 1 +V b =10m 3 -> Vi = (10-V b ) 

(10 -v b \p x+ v bPb 


Pl= 


(lO-F 4 )+r 4 


K = 0.6226 


>F,= 9.3774 m 3 


It is shown that 0.6226 m 3 of Barite requires 9.3774 m 3 of water to produce 10 m 3 mud of 
density 1.393 kg/1. 
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5.7 ECD. Barite 

Start by converting to SI units. When p mud > ~ 1.8 kg/1, the barite concentration is high. At stillstand, 
most mud types quickly develop a yield point, and barite is suspended in the gelled mud. During 
pumping of mud in laminar mode, any existing gel is broken, and barite will not be held in suspension as 
during stillstand with gelled mud. However, the mud will under these conditions exhibit a high effective 
viscosity but barite will, nevertheless, slip slowly downwards due to gravity. Agglomeration may also 
affect the settling rate. Barite has a specific density of 4.2, and an average particle size of approximately 
20 [im. In accordance with Stokes’s law of settling, barite particles will settle, but as indicated above, 
very slowly. However, wellbore inclination will shorten the settling distance to only a few centimetres, 
and a stratified bed of barite and cuttings will form. The stratified layers of solids will, when the angle 
of repose is surpassed, slide downwards. The critical angle of repose in found mainly in the build-up 
section of the well, and here sliding takes place. After the first “landslide”, the barite concentration in the 
mud has become lower and accordingly, also the viscosity. Barite settling will now be a fraction faster, 
and the process continues at accelerated speed. The self-preserving dynamics may lead to separation of 
almost all the barite in the fluid in the build-up section, especially during still stand after tripping-out 
of the hole. The pile of Barite at the bottom of the build-up section may cause pipe sticking. Worse is 
the reduction of mud weight. It may disrupt the safety. 

The following summary of guidelines to obtain efficient hole-cleaning is based on field experiences and 
on many laboratory investigations: 


a) Apply micro barite or equivalent. Such material systems have an average particle size of 2 [im 

b) Use of top-drive to allow for pipe rotation and redistribution of the settling barite and cuttings 
while tripping and drilling 

c) Maximize fluid velocity during drilling and reaming by increasing pumping output and/or use 
large diameter drill string. This improves hydraulic and agitation of cuttings. Use in addition 
downhole flow enhancers 

d) Design mud rheology so that it enhances turbulence in the annulus (if ECD allows). Turbulent 
flow will bring the cuttings and barite back into the flow stream if point 3 above “failed” 

e) Perform frequent wiper trips 

5.8 ECD. Flow rate & fluid consistency 

a) In order to determine the ECD we need to estimate the pressure loss in the annulus. The mud 
has some gel, so we assume the Bingham model will suit well. 




T 600 r 300 

Y 600 — ^300 


51.7-30.6 

1022-511 


= 0.0413 Pas 


t 0 =51.7 Pa -1022 • 0.0413 Pas = 9.5 Pa 
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Although the flow rate is rather high, we assume laminar flow in the annulus. However, if the 
purpose were to find the exact Ap we would need to check the Reynolds number to confirm 
this. Annular pressure loss: 


A P a = 48 */V 


Lpc/OH * V DC/OH + Lpp/OH ’ V DP/OH _j_ ^DP/C * V DP/C 

(d o — d t ) DC/OH (d Q — d t ) DP/OH ( d Q — d t ) DP/C 


+ 


6 -V 


DC/OH 


DP/OH 


^DP/C 


(d 0 -d i ) L 


(d 0 di) L 


(d„-d,) L 


= A p a = 385179 Pa - 3.85 bar 


Now the ECD: 


ECD = Phydr + A /V„„ = p- g’h + 
g-h g-h 

3 85 10 5 

ECD = 1250 +- : — : -= 1270 kg/m 3 

0.0981-2100 

b) What is the ECD during drilling? Production of mud and cuttings: 
q mud = 3 500 /(1000 • 60) = 0.058 m 3 /s 

Stings = R0P ' A bi,= 3 500/(1000 • 60) • 3.14/4 • (17.5 -0.0254) 2 = 0.0022 m 3 /s 
Total flow rate in annulus: q q , q ~ q , 

T-ann = T-mud + ^-cuttings = ^-mud 

Volumetric concentration of cuttings: 

c = q /(q ,) = 0.0022 / 0.058 = 0.038 

cuttings ^-cuttings v T-mud 7 

We assume slip; R = 0.75 

c = 0.038 / 0.75 = 0.051 

cuttings 

ECD contribution will be enhanced by the higher density of cuttings: 

ECD — p mu d ’ C mu( i + Pcuttings ’ C cuttings ~ 1-25(1 — 0.051) + 2.4 ■ 0.051 — 1.31 kg/l 


The effect of longer well is neglected. 

The casing shoe equivalent fracture density: p fr = p/ gh = 228.10 5 / (1 800.9.81) = 1 291 kg/m 3 
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Exercises in Drilling Fluid Engineering 


Wellbore challenges 


c) Suggested procedure to determine pressure profile, which is a f(T): 


1) Find the stable temperature profile in the well in accordance with Chapter 7 of the 
text book. Default temperature gradient in sedimentary rocks = 10 + 0.03°C/m. 

2) Look up the thermal coefficients of materials, a, in a handbook; 

OCsteel — 

OCwater 

3) p, = p a — p a ■«( ,■ T ) at each depth interval Ah 

4) Estimate Pi = ))) 

5) Find p bottom as J^P, 8 Al > 

s-\ / , _ Pbottom 

b) P eguiv - 

5.9 Water activity 

a) One definition of water activity is the vapor pressure of its salt water solution and compares 
it to the water vapor pressure of distilled water. 


A 


w 


P saline 


water 


P destiilled 


water 



qaiteye 

Challenge the way we run 


EXPERIENCE THE POWER OF 
FULL ENGAGEMENT... 


RUN FASTER. — r 

RUN LONGER.. 

RUN EASIER... > 
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Wellbore challenges 


Water solution: Distilled water thus has a water activity of 1.00. All water molecules are 
freely moving around, not electrostatically bound to salt ions. A w will diminish at increasing 
salt content. 

Pore water: Find the A w of clay through its in situ weight and then by comparing with clays 
with known A w . Thus A w of water phase is determined, and type of salt. 


b) Salt ions attract polar water and make it in-active. Distilled water has a water activity 1.0 
whereas for saline water it is less than 1.0. 

c) Clay swelling leads to increased pressure according to this equation: 


RT 

P swelling T/ r ^ 


(p j ^ 

Jr' w.clay 

V P w.mud J 


V w = Molar volume of water vapor; 


Most shale formations contain water since shale is (water wet). Water sensitive clay materials, 
such as smectite, illite and mixed - layer clays, will adsorb water and swell. This leads to an 
elevated localized pressure around the wellbore. The pressure may surpass the material strength; 
the material will disintegrate and collapse. Water is transported into the shale through different 
mechanisms (hydraulic pressure, osmosis, diffusion etc.). 


d) Adding K-salt to WBM. Dissolved salt will bind water and hinder it to evaporate; the water vapor 
pressure will declines and so also the measured A w . The A w is inversely proportional to salinity. 

e) If A /f ^ A , will be drawn from the OBM and into the formation, making the shale 
swell brittle. 

If A / , '\'d a a ^ A , the shale will become brittle. 

w, OBM w, pore 

Yes, it can be reduced by creating a filter in the shale. This can be done by adding particles to 
the water phase where the average particle size distribution — 1/3 of the average pore throat 
size distribution. 


5.10 Shale stability 

a) Using Figure 5-6 in the task section a calcium chloride concentration of 28.2% by weight is 
needed to give an activity of 0.69. 

b) Take an “in-situ” cut of the core, crush it and place it in a desiccator (completely air tight), 
expose it to different humilities (vapour pressure) through designed brines. When weight of the 
exposed sample has been converted to typical weight evolution of this type of shale. Through 
this comparison the water activity in the core can be determined. Since this is the water activity 
of the shales pore water, the water activity on the mud should be identical to avoid swelling 
or shrinking of the shale. 
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Wellbore challenges 


c) Oil does not penetrate clay because the capillary forces are too large; 
2• a f - cos# 


This depends on the wettability (which governs the 0). For oil to be the wetting phase against 
the formation/cuttings/steel a wetting agent must be added to the base oil. If water is emulsified 
into the oil, the water activity of the water phase must be controlled properly 

Emulsifier is used to reduce the surface tension between oil and water, enabling smaller droplets, 
enhancing stability Wetting agents make sure that oil is wetting the shale. 

1 Shale stability 

a) By adding salt to the muds water phase corresponding to A w of pore water. No osmotic force 
between mud and clay, and thus no water driven swelling. There are 5 ingredients in OBM, 
the A w and the salt type should be as close to pore waters as possible. 

1. Base oil 

2. Water 

3. Emulsifier (surface active additive I) 

4. Surface wetting (surface active additive II) 

5. Salt 

a) K + is geometrically suitable in between Montmorilonite platelets and leads to high platelet 
attraction (low swelling). 

b) The higher concentration, the more Na ++ is exchanged. 

c) Yes, it can be reduced by creating a filter in the shale. This can be done by adding particles to 
the water phase where the average particle size distribution — 1/3 of the average pore throat 
size distribution. 

d) One emulsifies, the other one oil-wets solid surfaces 

e) Through Cheneverts dissicator method 

f) Water flow, ion-flux, pressure flux 

g) Osmotic, non-invading drilling fluid 

h) Cyclic spalling, but also fractures and kick 
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Wellbore challenges 


5.12 Wellbore problems 

a) Lots of material (cuttings, cavings, sloughing shale) combined with small and large washouts 
where material tends to accumulate. During tripping, material is scraped and squeezed. 



Figure 5-7. Differential sticking mechanism. 

Differential sticking: Once the drill string touches the wellbore wall, the mud pressure, which is 
higher than the formation fluid pressure, holds the string in place. As mud cake builds around 
the string, a pressure seal is formed and sticking occurs as shown in Figure 5-7. 
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Wellbore challenges 


b) The consequences of a stuck pipe are very costly. They include: 

• Lost drilling time spent on freeing the pipe, and on fishing if not ask to jar out string. 

• Abandon tools in the hole because fishing was given up. 

Rules of Thumb 

1. Begin working the string immediately. Jar in the opposite direction to the pipe movement 
prior to becoming stuck. 

2. Work the pipe to the limits. 

3. If getting movement down, concentrate on expanding downward movement - and 
vice versa. 

4. In cases of hole bridging/packing off, concentrate on downward working. Increase 
applied working force in gradual increments up to the maximum. 

c) Salt tends to “flow” at high temperature and pressure. Use high mud weight and a mud that 
tend to dissolve salt (to keep the wellbore diameter as large as possible). 

WBM can be used with some additives. Halite has little creep tendency. 

d) Dominating mechanisms at lost circulation are: 

- low pressure window 

- high ECD (for many reasons) 

Counter measures are: 

- improved hole cleaning 

- improved ECD - control 
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Wellbore challenges 


Wellbore breathing (ballooning) (any two sentence is enough for get full marks) 

• The onset of wellbore breathing, often referred to as wellbore ballooning, is typically an indicator 
of imminent lost circulation. 

• Wellbore breathing is associated with fractures that open when annular pressure is applied to 
the Wellbore and close when the pressure is reduced. 

• These fractures fill with drilling fluid when open and subsequently return the fluid is observed 
as a flow out of the Wellbore when the pumps are off. 

• One of the more severe consequences of wellbore breathing is the misinterpretation of the 
observed flow as a kick when the pumps are shut down. 

• Implementing well-control procedures and increasing the mud weight is often enough to 
propagate the existing fractures leading to severe loss of circulation, a situation that is much 
more difficult to manage. 
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